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ABSTRACT 


The  parameters  of  Low  Probability  of  Intercept  (LPI)  radar  signals  are  hard  to 
identify  by  using  traditional  periodogram  signal  processing  techniques.  Using  the 
Wigner  Distribution  (WD),  this  thesis  examines  eight  types  of  LPI  radar  signals.  Signal 
to  noise  ratios  of  0  dB  and  -6dB  are  also  investigated. 

The  eight  types  LPI  radar  signals  examined  include  Frequency  Modulation 
Continuous  Wave  (FMCW),  Frank  code,  PI  code,  P2  code,  P3  code,  P4  code,  COSTAS 
frequency  hopping  and  Phase  Shift  Keying/Frequency  Shift  Keying  (PSK/FSK)  signals. 
Binary  Phase  Shift  Keying  (BPSK)  signals  although  not  used  in  modem  LPI  radars  are 
also  examined  to  further  illustrate  the  principal  characteristics  of  the  WD. 
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I.  INTRODUCTION 


A.  LPI  RADARS 

On  the  modem  battlefield,  radar  systems  face  serious  threats  from  EA  (Electronic 
Attack)  and  ARMs  (Anti  Radiation  Missiles).  To  perform  the  necessary  target  detection 
and  tracking  and  simultaneously  hide  themselves  from  an  enemy  attack,  radar  systems 
should  be  Low  Probability  of  Intercept  (LPI).  [Ref.  1] 

Low  probability  of  intercept  radars  have  the  property  of  low  power,  wide 
bandwidth,  frequency  variability  or  other  design  features  to  make  them  difficult  to  be 
detected  or  identified  by  passive  intercept  receiver  devices  such  as  Electronic  Support 
(ES)  or  Electronic  Intelligence  (ELINT)  receivers.  [Ref.  2]  LPI  radars  attempt  to  provide 
detection  of  targets  at  longer  range  than  intercept  receivers  can  accomplish  detection  of 
the  radar.  The  success  of  LPI  radars  is  measured  by  how  hard  it  is  for  them  to  be 
detected. 

The  LPI  radars  typically  have  the  following  features: 

(1)  Very  low  peak  power  (on  the  order  of  1W) 

(2)  Radiated  energy  spread  over  a  wide  angular  region,  over  a  long  time  interval, 
and  over  a  wide  frequency  band. 

(3)  Continuous  Wave  (CW)  radiation  with  a  large  time  bandwidth  product. 

(4)  Low  sidelobe  transmit  antenna  and  carefully  selected  carrier  frequency. 

(5)  Reduced  receiver  noise  temperature  and  loss.  [Ref.  1,  3,  4] 
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The  transmitter  makes  use  of  sophisticated  frequency  and  phase  modulation  to 
spread  the  signal  bandwidth  making  the  signal  hard  to  intercept.  The  receiver  makes  use 
of  the  appropriate  matched  filter  so  that  the  radar  perfonnance  is  similar  to  that  of 
traditional  pulsed  radar  radiating  the  same  amount  of  average  power. 

In  this  thesis,  Binary  Phase  Shift  Keying  (BPSK),  Frequency  Modulation 
Continuous  Wave  (FMCW),  polyphase  code  modulation  (Frank,  PI,  P2,  P3,  P4), 
COSTAS  frequency  hopping,  and  PSK/FSK  signals  will  be  analyzed  using  the  Wigner 
distribution  in  order  to  examine  detection  capability  and  to  extract  the  modulation 
parameters. 

B.  TIME-FREQUENCY  ANALYSIS  (TFA)  DETECTION  TECHNIQUES 

In  many  operation  radar  environments,  the  non- stationary  nature  of  the  received 
radar  signal  mandates  the  use  of  some  form  of  time-frequency  analysis  (TFA),  which  can 
clearly  bring  out  the  non-stationary  behavior  of  the  signal.  The  important  virtue  of  TFA 
is  that  it  provides  an  indication  of  the  specific  times  during  which  certain  spectral 
components  of  the  signal  are  observed.  [Ref.  5] 

Time-frequency  analysis  allows  time  and  frequency  resolution  to  be  controlled 
independently,  but  when  used  for  analysis  of  multi-frequency  component  signals  this 
approach  is  vulnerable  to  cross-terms  arising  midway  between  frequency  components. 
There  are  five  techniques,  which  are  in  current  use.  [Ref.  6] 

1.  Short-Time  Fourier  Transform  (STFT)  for  a  discrete  data  x(l )  is  defined  as 

oo  _2  •  I  n 

STFT {l, Jfc)  =  £  w*  (/  - n) ' 40 ' e*J  N  ( 1  •  1 ) 
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Estimating  of  the  TFA  using  the  STFT  involves  computing  a  Fourier  transform  (FT)  of  a 
discrete  time  series  x{l)  using  a  sliding  window  wN .  At  each  time  index  /,  a  N  point  FT 

is  computed,  giving  a  spectral  estimate  at  the  frequency  index  k .  The  window  length  N 

is  chosen  to  optimize  the  time  and  frequency  resolution;  the  signal  is  assumed  to  be 
stationary  within  the  window.  Signal  energy  density  can  be  estimated  from  the  squared 
modulus  of  the  STFT.  This  is  called  the  spectrogram.  [Ref.  6] 

SP(l,k)  =  \STFT(I,k)\2  (1.2) 

2.  Wigner  Distribution  of  a  continuous  signal  x(t)  is  defined  as 

co  __ 

WD(t,co)=  [  x(t  +  —  )-x*(t~—  )-e~j<aTdr  (1.3) 

J  2  2 

T=~  00 

where  <x>  is  angular  frequency  2nf ,  the  integral  is  from  -oo  to  oo ,  and  the  *  indicates  the 
complex  conjugate  of  the  signal  x(t)  .  Various  formulations  of  the  WD  are  possible  for 
sampled  time  series.  One  of  these  uses  a  moving  window  wN  with  2A-1  nonzero  values. 
[Ref.  6] 

3 .  Wigner  Distribution  (WD)  for  a  discrete  data  x(/)  is  defined  as 

V-l  -Ixjk— 

WD(l,k )  =  2  y,  wN(n)-x(l  +  n)-x*(l-n)-e  N  (1.4) 

n=—N 

The  WD  offers  very  good  time  and  frequency  resolution  compared  with  the  SP,  but  the 
cross-terms  generated  from  analysis  of  multi-component  signals  can  make  visual 
interpretation  of  the  TFA  difficult.  If  another  window  function  wM  with  2 M  - 1  nonzero 
values  is  introduced,  the  cross-terms  can  be  reduced.  [Ref.  6]  This  gives  the  smoothed 
WD. 


3 


4.  Smoothed  Wigner  Distribution  (SWD)  for  a  discrete  data  x(l)  is  defined  as: 

r  "-4  *  ]  -iKjk± 

SWD(l,k)  =  2  ^  wN(ri) \  ^  wM(m)-x(I +  m  +  n)-x  (I +  m-n)>-e  N  (1.5) 

n=—N  \<m=—M  J 

The  discrete  time  of  WD  (1.4)  is  periodic  in  n ,  but  for  SWD  the  discrete  time  is  periodic 
in  2 n  .  So  unless  the  signal  x(l)  has  been  sampled  at  least  twice  the  Nyquist  rate,  there 
will  be  aliasing  in  the  WD.  This  problem  can  be  solved  either  by  sampling  the  signal  at 
twice  the  Nyquist  rate  or  by  estimating  the  TFA  of  the  real  valued  time  series  x(I)  .  [Ref. 
6] 

5.  Choi- Williams  Distribution  (CWD)  for  a  discrete  data  x(l)  is  defined  as: 

AM  I  M- 1 

CWD(l,k)  =  2YjwN{n)\  X  wJm)- 

n=—N  m=—M 

(1.6) 

where  the  parameter  cr  can  be  adjusted  to  attenuate  cross-terms.  The  CWD  is  a  further 
refinement  of  the  SWD,  in  which  extra  smoothing  is  achieved  with  an  exponential 
weighting  function.  Lower  values  of  cr  attenuate  cross-terms  but  degrade  time  and 
frequency  resolution.  For  large  values  of  cr ,  cross-terms  are  not  suppressed  at  all,  and 
the  CWD  becomes  equivalent  to  the  SWD.  A  good  choice  of  cr  is  in  the  range  0.1  to  10. 
[Ref.  6] 


In  the  following  chapters,  the  WD  and  its  results  will  be  introduced  and  examined. 
The  SWD  and  the  CWD  are  for  the  further  studies  of  the  TFA  techniques.  They  will  not 
be  examined  in  this  thesis. 
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c. 


PRINCIPAL  CONTRIBUTIONS 


The  objective  of  this  thesis  was  to  develop  the  Wigner  distribution  (in 
MATLAB®  [Ref.  7])  and  use  it  to  analyze  and  extract  the  parameters  from  a  variety  of 

complex  LPI  radar  modulations.  In  order  to  use  the  Fast  Fourier  Transform  (FFT),  a 
special  kernel  transformation  is  developed  and  two  simple  examples  are  worked  in  order 
to  insure  that  the  kernel  transformation  is  calculated  correctly. 

The  114  signals  analyzed  include  BPSK,  FMCW,  Frank  code,  PI  code,  P2  code 
P3  code,  P4  code,  COSTAS,  PSK/FSK  with  frequency  hopping  and  two  test  signals. 
Using  the  Wigner  distribution,  this  thesis  successfully  extracted  all  the  parameters  of  the 
LPI  signals  except  the  PSK/FSK  with  target  frequency  hopping. 
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D.  THESIS  OUTLINE 

Chapter  II  presents  a  complete  introduction  of  the  Wigner  Distribution  (WD). 
Two  simple  examples  are  provided  in  this  chapter. 

Chapter  III  presents  the  Binary  Phase  Shift  Keying  (BPSK)  signals  and  the 
analysis  of  BPSK  signals  using  the  Wigner  distribution. 

Chapter  IV  presents  the  FMCW  signals  and  the  analysis  of  FMCW  signals  using 
the  Wigner  distribution. 

Chapter  V  presents  the  polyphase  signals  (Frank  code,  PI,  P2,  P3  and  P4)  and  the 
analysis  using  the  Wigner  distribution.  This  chapter  also  presents  the  comparison  of 
polyphase  signals. 

Chapter  VI  presents  the  COSTAS  frequency  hopping  signals  and  the  analysis  of 
COSTAS  signals  using  the  Wigner  distribution. 

Chapter  VII  presents  the  Phase  Shift  Keying/Frequency  Shift  Keying  (PSK/FSK) 
combined  with  frequency  hopping  signals  and  the  analysis  of  PSK/FSK  signals  using  the 
Wigner  distribution. 

Chapter  VIII  presents  some  concluding  remarks  and  recommendations. 

The  Appendix  contains  the  MATLAB®  m-fdes  of  the  Wigner  distribution  and  the 
LPI  radar  signal  matrix. 
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II.  THE  WIGNER  DISTRIBUTION 


The  Wigner  Distribution  (WD)  introduced  by  Wigner  in  1932  as  a  phase 
representation  in  quantum  mechanics  gives  a  simultaneous  representation  of  a  signal  in 
space  and  spatial-frequency  variables.  The  Wigner  distribution  has  been  noted  as  one  of 
the  more  useful  time-frequency  analysis  (TFA)  techniques  for  signal  processing. 

A.  WIGNER  DISTRIBUTION  EQUATION 

Recall  (1.3),  the  Wigner  distribution  of  a  continuous  input  signal  x(t)  is  defined 
as 

oo 

W{t,co)  =  |  x(t +  t  1 2)- x*  (t-r  1 2)e~i<oz  dr  (2.1) 

-oo 

where  t  is  the  time  variable  and  a>  is  the  frequency  variable.  The  Wigner  distribution  is 
a  two-dimension  function  describing  the  frequency  content  of  a  signal  as  a  function  of 
time. 

This  continuous  time  and  frequency  representation  can  be  modified  for  the 
discrete  sequence  x{l ) ,  where  /  is  a  discrete  time  index,  /  =  ...,-1,0,  1,  . .  .The  discrete 
Wigner  distribution  (WD)  is  defined  as 

oo 

W (I, o?)  =  2  ^  x(l  +  n)x\l -  n)e-j2,J>n  (2.2) 

«=— oo 

Further  modification  results  in  the  pseudo-Wigner  distribution  or  windowed-Wigner 
distribution,  which  is  defined  in  (2.3),  the  same  as  (1.4). 

N- 1 

W(l, of)  =  2  ^  x(l  +  ri)x* (/  - n)w(n)w(-n)e~J2a>n  (2.3) 

n=—N 
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where  x(l )  is  a  discrete  input  signal  with  /  from  -oo  to  oo ,  w(n)  is  a  length  2N real 
window  function  with  w(0)  =  1 .  Here  N  must  be  as  large  as  possible  within  the  limits  of 
an  acceptable  computational  cost  because  a  large  N  gives  more  output  samples,  yielding  a 
smoother  result.  Considering  the  rectangular  window  function  with  a  magnitude  equal  to 
one  the  WD  becomes  [Ref.  8] 

JV— 1 

W(l, co)  =  2  Y,  *(/  +  n)x\l  - n)ejla,n  (2.4) 

n=—N 

Using  ft{n)  to  represent  the  kernel  function 

f(n)  =  x(l  +  n)x*(l-n)  (2.5) 

the  WD  becomes 


JV— 1 

W (l,  a?)  =  2  ^  f(n)e~J2' 

n=—N 


where  the  continuous  frequency  variable  co  is  sampled  by  [Ref.  9,  10] 


(2.6) 


(2.7) 


and  where  k  =  0,1, 2, ...2 N  - 1 .  The  kernel  indexes  are  modified  to  fit  the  standard 
Discrete  Fourier  Transform  (DFT).  Since 

fl(n)  =  f*(~n)  (2.8) 


the  kernel  is  a  symmetric  function,  so  the  Discrete  Fourier  Transform  (DFT)  of  the  kernel 
is  always  real. 

From  equation  (2.6)  and  (2.7),  the  WD  becomes 
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(2.9) 


2  N  2  N 


Adjusting  the  limits  of  n  in  order  to  use  the  standard  FFT  algorithms,  (2.9)  becomes 


W(l,—)  =  2  Y  fi(n)  exp(-^^-) 
2  N  F  2  N 


(2-10) 


Note  that  in  (2.10)  the  kernel  function  has  been  adjusted  to  ft  (n) ,  where 


fi  (»)  =  i 


UM-2N) 


0  <  n  <  N  - 1 
n  =  N 


N  +  \<n< 2N 


(2.11) 


The  resulting  WD  is 


2iV-l 

W(l,co)  =  2YJMn)e~i2wn 

n= 0 


(2.12) 


2 

W(l,k)  =  2^f  (n) exp (-j  ——) 

to  N 


(2.13) 


Equation  (2.13)  is  the  final  WD  equation,  and  will  be  used  on  the  following  paragraph  to 
show  the  computational  procedure  of  the  WD  values. 


B.  TWO  SIMPLE  EXAMPLES 

1.  Example  With  Real  Inputs 

Suppose  there  is  an  input  signal 


x(/)  =  {2,4, 3, 6, 1,7} 

=  {2(1  =  -3),  4(1  =  -2),  3(1  =  -1),  6(1  =  0),  1  (/  =  1),  1(1  =  2)} 


(2.14) 


9 


where  the  length  of  the  input  signal  x(l)  is  2N  =  6,  or  N=3.  /  is  a  discrete  time  index 


from  -N toN-l.  Note  that  x  =  0  for  /  <  -4  or  / >  3  .  Recall  (2.1 1),  since  now  N=3,  so 


f,(n)  ,  0  <n  <2 
0  ,  n  =  3 


[fin- 6),  4  <  «  <  5 


(2.15) 


For  example,  for  1=  - 3,  the  kernel  f_3  (n)  becomes 

f-iin)  =  {f_3(n  =  0),f_3(n  =  1  ),f_3(n  =  2 ),f_3(n  =  3 ),f_3(n  =  4 ),f_3(n  =  5)} ,  from  (2.15) 


f-3  (n  =  3)  =  /_3 (n  =  3)  =  0 ,  /_3'(n  =  4)  =  f_3(n  =  4-6)  =  f_3(n  =  -2)  and 


f~3  (n  =  5)  =  f_3  (n  =  5  -  6)  =  /_3  (n  =  -1) 

Hence, 

/-3’(«)  =  {f-:M  =  0  ),/_3(«  =  l),/_3(n  =  2),0,/3(/z  =  -2),/_3(«  =  -1)} 


Recall  that  /;  («)  =  x(/  +  n)  -x*  (l-  n)  , 


and  x(/)  =  {2(/  =  -3),  4(/  =  -2),  3(/  =  - 1),  6(/  =  0),  1(/  =  1),  7(7  =  2)} 
and  x(/)  =  0  for  /  <  -4  or  /  >  3  . 

/_3(«)  for  input  signal  x(/)is  computed  as  follows 


f_3  ( n  =  0)  =  x(-3  +  0)  •  jc*  (-3  -  0)  =  x(-3)  •  x*  (-3)  =  2-2  =  4 
f_3  (n  =  1)  =  jc(-3  + 1)  •  x*  (-3  - 1)  =  x(-2)  •  x*  (-4)  =  4-0  =  0 
f_3  (n  =  2)  =  x(-3  +  2)  •  x*  (-3  -  2)  =  x(- 1)  •  x*  (-5)  =  3-0  =  0 
f~iin  =  3)  =  0 

/  3  («  =  4)  =  /_3  (n  =  -2)  =  x(-3  -  2)  •  x*  (-3  +  2)  =  x(-5)  •  x*  (- 1)  =  0  •  3  =  0 
/  3  (n  =  5)  =  /_3  (//  =  -1)  =  x(-3  - 1)  •  x*(-3  + 1)  =  x(-4)  •  x*(-2)  =  0-4  =  0 
So,  /_3'(/i)  =  {4, 0,0, 0,0,0} 
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And  similarly  for  /  =  0, 


f0  (n  =  0)  =  jc(0  +  0)  ■  x*  (0  -  0)  =  jc(0)  ■  /(O)  =  6  ■  6  =  36 
fQ  (n  =  1)  =  x(0  + 1)  ■  x*  (0  - 1)  =  jc(1)  •  jc*(-1)  =  1  •  3  =  3 
/„  (n  =  2)  =  jc(0  +  2)  •  /  (0  -  2)  =  x(2)  ■  x  (-2)  =  7-4  =  28 
fo(n  =  3)  =  0 

fo  in  =  4)  =  f0  (n  =  -2)  =  x(0  -  2)  •  (0  +  2)  =  x(-2)  ■  x  (2)  =  4  ■  7  =  28 

f0(n  =  5)  =  f0(n  =  -1)  =  jc(0  - 1)  •  x*(0  + 1)  =  jc(-1)  -  jc*(1)  =  3-1  =  3 
So,  f0(n)  =  { 36,3,28,0,28,3} 

Repeating  the  above  procedures,  the  kernel  matrix  for  /  =  -4  to  3,  and  n  =  0  to  5  is  as 
shown  in  Figure  1. 
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Figure  1 .  The  Kernel  f]  (n)  Matrix  for  the  Real  Six  Input  Example. 
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The  second  step  is  to  use  (2.13)  to  calculate  the  Wigner  distribution.  As  an  example  of 
the  calculation,  one  can  pick  any  /  and  k  to  examine  the  values  inside  the  WD  matrix. 

For  example,  choose  /  =  1,  k  =  2  and  N=  3.  The  WD  is 

,  Trial 

W(l=l,k=2)  =  2'£Jf,  (n)exp(-/— ) 

to  N 

2-3-1  ,  -2,7 

= 2  Z  /i  (w)  cxp(./  — i-)  (2- 1 6) 

w=0  ^ 

5  i  0  77"  M 

=  2Z/i(»)exp(-y— ) 

n=0  3 

From  the  kernel  matrix  in  Figure  1,  the  kernel  function  for  /  =  1  is 

f  (n)  =  {1,42,0,0,0,42} .  Which  means  that 

/'( 0)  =  1,//(1)  =  42,/;  (2)  =  0,/  (3)  =  0,/  (4)  =  0,/  (5)  =  42  . 

So  from  (2. 16),  the  WD  for  /  =  1,  k  =  2  is 

W  (■• 2)  =  ^ ' eXP<-2  > +  W‘> '  exp(-yT/i)  +  2/(2)  ■  exp(-,/pi) 

+  2,/;'(3)  •  exp(-,  +  2/;'  (4)  •  exp(-,'  +  2/(5)  •  exp,-, 

=  2 •  1  •  (0)  +  2 •  42 •  (-0.5000  -  0.8660i)  +  2-  0  +  2-  0  +  2-  0  +  2-42- (-0.5000  +  0.8660i) 
=  -82 

Repeating  the  above  procedures,  gives  the  WD  matrix  at  discrete  time  index  /  =  -4  to  3 
and  discrete  frequency  index  k  =  0  to  5  as  a  symmetric  matrix  shown  in  Figure  2.  An 
important  feature  in  this  Wigner  distribution  computation  result  is  that  all  the  components 
in  WD  matrix  are  real. 
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Figure  2.  The  WD  Matrix  W(l,k )  for  the  Real  Six  Input  Example. 

2.  Example  With  Complex  Inputs 

Consider  the  input  signal 

x  =  I  +  j-Q 


where 

/  =  cos(2;r  fct ) 
Q  =  sm(2nfct) 


(2.17) 


(2.18) 

(2.19) 


If  carrier  frequency  fc= 1kHz;  sampling  frequency  f  =  7kHz;  and  t  =  0 : 


— :  — ,  then  the 

fs  fs 


first  eight  input  points  from  discrete  time  index  /  =  -4  to  3  is 
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x(I)  =  {1  +  0/,  0.6235  +  0.7818/,  -0.2225  +  0.9749/,  -0.9010  +  0.4339/, 
-0.9010  -  0.4339/,  -0.2225  -  0.9749/,  0.6235  -  0.7818/,  1  +  0/} 


Consider  the  value  when  l  =  0,  n  =  3. 

Using  (2.11)  with  an  input  length  2N=  8  or  N  =  4.  The  kernel  is 


[fiin- 8) 


,  0  <  n  <  3 
,77  =  4 
,  5  <  77  <  7 


(2.21) 


Here  // (77)  =  {/  (1),  / (2), / (3), 0,  / (-3), / (-2),  f, (-1)} 

Since  ft  (77)  =  x(/  +  77)  -x*(l-  n)  ,  the  kernel  at  /  =  0,  n  =  3  is 

/o' (3)  =  /0(3)  =  x(0  +  3)  •  x  (0  -  3)  =  x(3)  •  x*(-3) 

=  1-  (0.6235  +  0.7818-/)* 

=  0.6235-0.7818  / 

Repeating  the  same  procedures  as  discussed  in  real  input  case,  the  kernel  matrix  for  the 
eight  inputs  complex  input  example  is  shown  in  Figure  3. 
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Figure  3.  The  Kernel  Matrix  for  the  Complex  Eight  Input  Example. 


Now  the  WD  matrix  is  examined  from  Figure  3,  when  /  =  -1,  the  kernel  is 

—  {1.0001,  -0.2225  +  0.9749/,  -0.9009 -0.4339/,  0.6235-0.7818/, 
0,0.6235  +  0.7818/,  -0.9009  +  0.4339/,  -0.2225-0.9749/} 


Consider  the  case  when  /=  -1,  k=  4.  Recall  (2.13),  the  WD  for  N  =  4  is 

2J V-l  , 

W(l  =  -\,k  =  4)  =  2YJfi  («) exp(-/  — -) 

to  N 

2£-i  ,  jrAn 

=  2i/_1(»)^)  (2-23) 

n= 0  4 

7 

=  2^  f-i(n)  exp(-jn^r ) 

n= 0 


From  (2.22)  and  (2.23) 
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w (1  =  -1,  k  =  4)  =  2  Yjf-i(n)  ■  exp irjrm) 

n= 0 


f- 1  (0) '  exp(0)  +  /’,  (1)  •  exp(-j^)  +  /',  (2)  •  (-/2^)  +  /',  (3)  •  (-j3;r) 
(4) •  (-/4/r)  +  /’  (5) •  (~j57r )  +  /’ (6) •  (-j6;r)  +  /',  (7) •  (-/ 7^) 


=  2- 


1 .000 1  •  (1)  +  (-0.2225  +  0.97490  ■(-!)  +  (-0.9009  -  0.43390  •  (1) 
+(0.6235  -  0.78 1 80  •  (-1)  +  0  +  (0.6235  +  0.78 1 8)  •  (-1) 
+(-0.9009  +  0.4339/)  ■  (1)  +  (-0.2225  -  0.9749/)  •  (-1) 


=  -3.2073 

Again  the  WD  matrix  of  the  eight  input  points  is  another  real  and  symmetric  2Nx2N 
matrix.  Note  the  important  feature:  The  WD  is  always  real  whether  the  input  signal  is 
real  or  complex.  Figure  4  shows  the  WD  matrix  of  the  eight  complex  input  samples. 
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Figure  4.  The  WD  Matrix  for  the  Complex  Eight  Input  Example. 
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Figure  5(a)  shows  the  mesh  plot  of  the  eight  input  complex  example  in  time 
domain.  The  time  resolution  1 /  fs  is  indicated  on  this  plot.  Figure  5(b)  shows  the  mesh 
plot  in  frequency  domain.  The  carrier  frequency  in  this  plot  is  about  900Hz,  very  close  to 
real  value  1kHz.  The  frequency  resolution  /'  12 1#  samples  is  also  indicated  on  this  plot. 


Time  Resolution  =  l/fs=  1/7000  =  0.1429  ms 

Mesh  of  8  inputs  complex  example 


1kHz  Carrier  frequency  Freq.  Resolution  =  fs/2/#samples 


(b) 

Figure  5.  WD  for  Eight  Inputs  Complex  Example  (a)  2D  Mesh  in  Time  Domain,  (b)  2D 

Mesh  in  Frequency  Domain. 
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Figure  6(a)  shows  the  3D  mesh  plot  of  the  eight  inputs  complex  example.  This 
plot  shows  the  magnitude  in  both  time  and  frequency  domain.  Figure  6(b)  shows  the 
contour  plot.  The  contour  plot  is  a  2D  time-  frequency  domain  plot.  The  magnitude  is 
represented  by  different  color  as  shows  in  the  color  bar. 


Mesh  of8  inputs  complex  example 


(b) 

Figure  6.  WD  for  Eight  Inputs  Complex  Example  (a)  3D  Mesh,  (b)  Contour. 
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Figure  7(a)(b)  shows  the  2D  time  domain  and  frequency  domain  mesh  plot  of  the 


101  inputs  complex  example  (Eq.  (2. 18),  (2.19)  with  t  =  0 :  — : - ).  The  1kHz  carrier 

•f S  f s 

frequency  is  more  clearly  shown  in  Figure  7(b)  than  that  in  Figure  5(b).  Comparing  with 
Figure  5  and  Figure  7  one  can  find  the  more  input  sample  points  that  are  used  the  better  is 
the  WD  output  performance. 


Mesh  of  101  inputs  complex  example 


(a) 


Figure  7. 


Mesh  of  101  inputs  complex  example 


Frequency  (Hz) 


(b) 

WD  for  101  Inputs  Complex  Example  (a)  2D  Mesh  in  Time  Domain,  (b)  2D 
Mesh  in  Frequency  Domain. 
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Figure  8(a)  shows  a  zoom  in  of  the  mesh  plot  in  time  domain.  The  total  time 
interval  and  time  resolution  1 1  fs  are  indicated  on  this  plot.  Figure  8(b)  shows  the  zoom 
in  mesh  plot  in  the  frequency  domain  and  indicates  the  frequency  interval  and  frequency 
resolution  fj  2  /#  samples  . 


Time  resolution  =  fs 

_  l 

7000 


Time  interval 

—  -i-*  (#points) 


Freq.  resolution  =  2 -^samples  Frefl  interval 

=  fs!  2=  35  00  Hz 


Figure  8.  WD  for  101  Inputs  Complex  Example  (a)  Zoom  In  2D  Mesh  in  Time  Domain. 

(b)  Zoom  In  2D  Mesh  in  Frequency  Domain. 
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Figure  9(a)  shows  the  three-dimension  mesh  plot  of  the  101  inputs  complex 
example.  This  plot  shows  the  magnitude  in  both  time  and  frequency  domain.  Figure  9(b) 
shows  the  contour  plot  of  the  101  inputs  complex  example.  Contour  plot  is  a  2D  time- 
frequency  domain  plot.  The  magnitude  is  represented  by  different  color  as  shows  in  the 
color  bar.  The  carrier  frequency  is  also  marked  on  the  plot. 


Mesh  of  101  inputs  complex  example 


-50 

0.015 


4000 


Frequency  (Hz) 


(a) 

Contour  of  101  inputs  complex  example 


(b) 

Figure  9.  WD  for  101  Inputs  Complex  Example  (a)  3D  Mesh,  (b)  Contour. 
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Now  consider  the  WD  for  two  carrier  frequencies  with  fcl  =  1  kHz ,  fc2  =  2 kHz  as 


shown  in  Figure  10.  Recall  that  /  =  cos(2;r  fct)  and  Q  =  sin(2;r  ft) .  Now 


/  =  cos(2;r fcXt )  +  cos(2;r fc2t )  and  Q  =  sin(2;r fclt)  +  sin(2;r fc2t ) 

Figure  10  (a)  is  the  2D  frequency  domain  mesh  plot  which  shows  the  two  carrier 
frequencies  as  well  as  the  frequency  cross  term.  The  frequency  of  the  cross  term  is 


fcl  +  fcl 


2 


.  Note  that  the  magnitudes1  are  positive  for  the  two  real  carrier  frequencies.  But 


for  the  cross  term,  the  magnitude  includes  both  the  positive  and  negative  parts. 


Figure  10(b)  shows  the  contour  plot  of  these  two  carrier  frequencies  and  their 
cross  term.  Note  that  the  shape  and  color  of  the  cross  term  is  not  like  those  of  the  two 
real  carrier  frequencies.  The  cross  term  contains  many  red  and  blue  ‘dots’,  but  the  real 
fd  fd  are  yellow  ‘strips’. 


Figure  11(a)  shows  the  3D  mesh  plot.  In  this  plot  displays  the  cross  term  is 
combined  by  many  ‘peaks’.  The  positive  peaks  represents  the  red  dots,  the  negative 
peaks  represents  the  blue  dots  in  figure  10(b). 

Figure  1 1(b)  is  the  2D  mesh  in  the  time  domain.  It  shows  the  real  signal  and  the 
cross  tenn  in  the  time  domain.  This  plot  states  that  the  cross  term  is  consisted  of  a  series 
of  positive  and  negative  magnitude  components  in  time  domain. 


1  The  label  “Magnitude”  was  used  for  WD  value  in  all  the  mesh  plots.  Note  that  is 
not  a  true  magnitude. 
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Mesh  of  two  carrier  frequencies  example 


(a) 


Contour  of  two  carrier  frequencies  example 


(b) 

Figure  10.  WD  for  Two  Carrier  Frequencies  Example  (a)  2D  Mesh  in  Frequency 

Domain,  (b)  Contour. 
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Mesh  of  two  carrier  frequencies  example 
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Mesh  of  two  earner  frequencies  example 


(b) 


Figure  11.  WD  for  Two  Carrier  Frequencies  Example  (a)  3D  Mesh,  (b)  2D  Mesh  in  Time 

Domain. 
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C.  SUMMARY 

This  chapter  contains  some  simple  example  signals  to  generate  the  kernel  function 
and  Wigner  distribution  (WD).  No  matter  whether  the  signals  are  real  or  complex,  the 
kernel  and  WD  matrix  are  always  real  and  symmetric.  This  is  a  very  important  feature 
for  the  Wigner  distribution  and  an  important  reason  why  the  WD  is  a  good  tool  to  analyze 
signals  by  the  mesh  and  contour  plots. 

In  Figures  5  through  11,  the  carrier  frequencies  are  very  easy  to  identify  in  these 
simple  examples.  In  the  following  chapters,  the  signals  contain  complex  LPI 
modulations  and  are  no  longer  that  simple.  The  Wigner  distribution’s  ability  to 
determine  the  phase  and  frequency  modulations  for  those  signals  will  be  demonstrated. 

Figures  10  and  11  illustrate  the  frequency  cross  terms  and  how  they  are 
represented  in  the  WD  mesh  and  contour  plots.  It  is  shown  that  the  cross  terms 
sometimes  make  the  signal  analysis  more  difficult. 

In  the  next  chapter,  the  BPSK  signals  and  its  WD  analysis  will  be  examined. 
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III.  BINARY  PHASE  SHIFT  KEYING  (BPSK) 


A.  BPSK 

Phase  shift  keying  (PSK)  is  an  angle-modulated,  constant-amplitude  digital 
modulation.  PSK  is  similar  to  conventional  phase  modulation  except  that  with  PSK  the 
input  signal  is  a  binary  digital  signal  and  a  limited  number  of  output  phases  are  possible. 

Binary  phase  shift  keying  (BPSK)  has  two  phases  for  a  single  carrier  frequency 
(“binary”  meaning  “2”).  One  output  phase  represents  logic  1  and  the  other  is  logic  0.  As 
the  input  digital  signal  changes  state,  the  phase  of  the  output  carrier  shifts  between  two 
angles  that  are  180°  out  of  phase.  BPSK  is  a  form  of  suppressed-carrier,  square- wave 
modulation  of  a  continuous  wave  (CW)  signal.  [Ref.  1 1] 

Figure  12  shows  the  BPSK  modulator  phase  diagram.  Figure  13  shows  the  BPSK 
modulator  output  phase  vs.  time  relationship.  Normally,  there  are  many  carrier  cycles  per 
modulation  period. 


cos  <x>ct 
(+90°) 

-5  in  coct 

(180°)  <  | 
Logic  0 

-cos  coct 
(-90°) 


5  in  cot 
*  (0°) 
Logic  1 


Figure  12.  BPSK  Modulator  Phase  Diagram. 
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Figure  13.  BPSK  Modulator  Output  Phase  vs.  Time  Relationship. 

1.  Binary  Phase-Coded  Pulse  Compression 


Phase  changes  can  be  used  to  increase  the  signal  bandwidth  of  a  long  pulse  for 
purpose  of  pulse  compression.  A  long  pulse  of  duration  T  is  divided  into  N  subpulses 
each  of  width  th .  An  increase  in  bandwidth  is  achieved  by  changing  the  phase  of  each 
subpulse  (since  the  rate  of  change  of  phase  with  time  is  a  frequency).  A  common  form  of 
phase  change  is  binary  phase  coding,  in  which  the  phase  of  each  subpulse  is  selected  to 
be  either  0  or  n  radians  according  to  some  specified  criterion. 

The  radar  receiver  correlates  the  received  pulses  with  the  known  transmitted 
pulse,  usually  with  a  matched  filter.  Thus  output  of  the  matched  filter  will  be  a 
compressed  pulse  of  width  th  and  will  have  peak  N  times  greater  than  that  of  the  long 

pulse.  The  pulse  compression  ratio  equals  to  the  number  of  subpulses  N  =  T  /  tb~  BT , 
where  the  bandwidth  B  ~  1  /  tb. 
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The  matched  filter  output  extends  for  a  time  T  on  either  side  of  the  peak  response. 
The  unwanted  but  unavoidable  portions  of  the  output  waveform  other  than  the 
compressed  pulse  are  known  as  time  sidelobes.  When  the  selection  of  the  phase  is  made 
at  random,  the  expected  maximum  (power)  sidelobe  is  about  2 IN  below  the  peak  of  the 
compressed  pulse.  [Ref.  4] 

2.  Barker  Codes 

One  family  of  binary  phase  codes  that  produce  compressed  waveform  with 
constant  sidelobe  levels  equal  to  unity  is  the  Barker  code.  Figure  14  illustrates  this 
concept  for  a  Barker  code  of  length  seven.  There  are  only  nine  known  Barker  codes  that 
share  this  unique  property.  They  are  listed  in  Table  1.  [Ref.  12] 


(b) 

Figure  14.  Barker  Code  of  Length  7  (a)  A  Long  Pulse  with  7  Equal  Subdivisions  Whose 
Individual  Phase  are  0°(+)  or  180°(-).  (b)  Autocorrelation  Function  for  Plot  (a). 
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In  general,  the  autocorrelation  function  (which  is  an  approximation  for  the 
matched  filter  output)  for  code  length  N,  Barker  code  will  be  2 Ntb  wide.  The  main  lobe 

is  2 tb  wide.  The  peak  value  is  equal  to  N.  There  are  (N- 1)/2  sidelobes  on  either  side  of 
the  main  lobe.  Note  that  in  Figure  14(b)  the  main  lobe  is  7  and  all  the  other  sidelobes  are 
one.  [Ref.  12] 


Code  Length 

Code  Elements 

Sidelobe  Level,  dB 

2 

+  -,  +  + 

-6.0 

3 

+  +  - 

-9.5 

4 

++-+,+++- 

-12.0 

5 

+  +  +  -  + 

-14.0 

7 

+++ — +- 

-16.9 

11 

+  +  + - +  —  +  - 

-20.8 

13 

+++++ — ++-+-+ 

-22.3 

Table  1.  Barker  Codes. 


The  most  sidelobe  reduction  offered  by  a  Barker  code  is  -22.3dB  (13  bit  code 
length),  which  may  not  be  sufficient  for  the  desired  radar  application.  However,  Barker 
codes  can  be  combined  to  generate  much  longer  codes.  In  the  following  example  and 
Figure  15,  Barker  code  B54  is  a  Barker  code  combines  code  length  5  and  code  length  4. 
But  unfortunately,  the  sidelobes  of  a  combined  Barker  code  autocorrelation  function  are 
no  longer  equal  to  unity. 

B54  =  {  11101,  11101,00010,  11101  }  (3.1) 
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Figure  15.  Combined  Barker  Code  with  Code  Length  5  and  4. 


Figure  16  below  shows  a  basic  block  diagram  of  the  transmitter  design.  The 
signal  x(t)  is  a  CW  sinusoid,  which  for  this  project  was  set  to 

x(t )  =  sin(2;r  fct)  (3 .2) 

with  carrier  frequency  fc .  The  discrete  signal  xnT  is  created  by  sampling  x(t)  with  at  a 
sample  rate  at  least  twice  that  of  the  highest  frequency  of  the  CW  signal  to  obey  the 
Nyquist  sampling  theorem  and  avoid  aliasing.  The  sampling  frequency,  fs ,  is  the 

product  of  the  carrier  frequency  fc  and  the  sample  ratio  variable 

SAR  =fjfc  (3.3) 

that  is  set  by  the  user.  To  obey  the  Nyquist  sampling  theorem,  SAR  must  be  greater  than 
two  to  avoid  anomalies  that  may  arise  from  floating  point  operations  right  at  the  Nyquist 
rate.  Due  to  these  limitations,  the  minimum  useful  sampling  frequency  is  three  times  the 
signal  frequency,  or  SAR  >=  3.  [Ref.  13] 

Although  there  are  only  13  bits  in  this  Barker  code,  the  actual  sequence  of  bits  in 

brkseq  also  depends  on  the  sampling  frequency  and  the  number  of  complete  periods  of 
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signal  xnT  that  are  to  occur  within  the  same  timeframe  as  one  Barker  code  bit.  If  the 
number  of  periods  per  Barker  bit  ( NPBB )  is  set  to  one,  one  full  period  of  the  carrier 
frequency  fits  within  one  bit  of  the  13-bit  Barker  code. 


Figure  16.  BPSK  Transmitter  Block  Diagram. 


In  this  chapter,  twelve  BPSK  signals  will  be  examined.  They  are  numbered  from 
one  to  twelve  and  listed  in  Table  2.  A  1kHz  carrier  frequency  and  7kHz  sampling 


frequency  are  used  for  these  signals. 


No. 

BPSK 

Barker  Code 
Length 

NPBB 

SNR 

1 

B  1  7  7  1  s 

7 

1 

Signal  Only 

2 

B  1  7  7  1  0 

7 

1 

OdB 

3 

B  1  7  7  1  -6 

7 

1 

-6  dB 

4 

B  1  7  11  1  s 

11 

1 

Signal  Only 

5 

B  1  7  11  1  0 

11 

1 

OdB 

6 

B  1  7  11  1  -6 

11 

1 

-6  dB 

7 

B  1  7  7  5  s 

7 

5 

Signal  Only 

8 

B  1  7  7  5  0 

7 

5 

OdB 

9 

B  1  7  7  5  -6 

7 

5 

-6  dB 

10 

B  1  7  11  5  s 

11 

5 

Signal  Only 

11 

B  1  7  11  5  0 

11 

5 

OdB 

12 

B  1  7  11  5  -6 

11 

5 

-6  dB 

Table  2.  BPSK  Signals. 
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B.  WIGNER  DISTRIBUTION  FOR  BPSK 

Consider  the  same  twelve  BPSK  signals  in  Table  2  and  use  them  as  the  inputs  for 
Wigner  distribution  (WD).  The  WD  results  of  these  twelve  signals  are  shown  in  Figure 
17  through  20.  The  mesh  plots  show  the  frequency  domain  of  the  BPSK  signals  after  the 
WD  processing.  The  contour  plots  show  both  the  frequency  domain  and  time  domain  of 
the  results. 

In  Figure  17,  the  carrier  frequency  fc  can  be  clearly  identified  by  the  location  of 

the  highest  or  lowest  peak  value  in  Figure  17(a),  or  by  the  center  of  the  symmetric  shapes 
in  Figure  17(b). 

Secondly,  the  peak  magnitude  in  Figure  175(a)  is  about  600.  So  the  3dB 
bandwidths  should  be  the  frequency  range  between  300  on  both  sides,  which  extends 
from  500Hz  to  1500  Hz  in  Figure  17(a).  The  3dB  bandwidth  is  also  the  value  of  fc 

divided  by  NPBB.  That  is  equal  to  1000Hz  in  Figure  17,  18  ( NPBB  =  1),  or  200Hz  in 
Figure  19,  20  ( NPBB  =  5).  The  3dB  bandwidth  is  also  called  the  chip  rate. 

Thirdly,  in  Figure  17(a),  if  we  look  closely  within  the  3dB  bandwidth,  one  can 
find  there  are  15  peaks  in  the  bandwidth  range.  In  other  words,  there  are  14  intervals  in 
the  range  from  500Hz  to  1500Hz,  and  14  is  two  times  the  Barker  code  length.  In  Figure 
18(a),  22  intervals  are  in  the  3dB  bandwidth  range  since  the  Barker  code  length  is  eleven. 

Using  the  same  analysis  for  Figures  17  to  20,  we  can  find  the  carrier  frequency, 
3dB  bandwidth  and  barker  code  length.  However,  for  NPBB  more  than  one,  the  analysis 
of  Barker  code  length  becomes  difficult. 
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c. 


SUMMARY 


The  results  of  Wigner  Distribution  analysis  for  BPSK  signals  are  discussed  for 
two  situations.  When  the  NPBB  (the  number  of  cycles  per  barker  bit)  is  equal  to  one,  the 
Barker  code  length  can  be  easily  found  by  the  Wigner  Distribution.  But  the  disadvantage 
is  the  carrier  frequency  information  becomes  more  difficult  to  determine  in  low  SNR  (the 
result  figures  are  not  contained  in  this  thesis).  Second,  when  the  NPBB  is  higher  than 
one,  the  carrier  frequency  and  3  dB  bandwidth  can  be  easily  found  in  low  SNR,  but  the 
Barker  code  length  becomes  hard  to  analyze. 

Table  3  below  shows  the  detection  effectiveness  of  the  carrier  frequency,  3dB 
bandwidth,  code  period  and  Barker  code  bit  from  Figure  17  to  20.  They  are  done  by  the 
visible  inspection  of  the  WD  results. 

In  the  next  chapter,  the  FMCW  signals  and  their  WD  analysis  will  be  examined. 
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BPSK 

NO. 

Carrier  Frequency 

Bandwidth 

Code  Period 

Bits/Code 

1 

100% 

100% 

100% 

100% 

2 

100% 

110% 

100% 

100% 

3 

0% 

0% 

0% 

0% 

4 

100% 

100% 

100% 

100% 

5 

100% 

90% 

100% 

100% 

6 

0% 

0% 

0% 

0% 

7 

100% 

100% 

0% 

0% 

8 

100% 

120% 

0% 

0% 

9 

100% 

150% 

0% 

0% 

10 

100% 

100% 

0% 

0% 

11 

100% 

120% 

0% 

0% 

12 

100% 

150% 

0% 

0% 

BPSK  effectiveness  on  detection 


150% 


130% 


110% 


90% 


<u 
o> 
n 
-*-< 

o  70% 

<u 

Q. 

50% 


30% 


10% 


□  Only  signal 
■  0  dB 

□  (-)  6  dB 


Carrier  freq. (Hz)  Bandwidth  (Hz)  Code  period  (ms) 

Parameters 


Bits/code 


Table  3.  WD  Detection  Effectiveness  for  BPSK  Signals. 
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IV.  FREQUENCY  MODULATION  CONTINUOUS  WAVE  (FMCW) 


A.  FMCW 

FMCW  radar  is  a  technique  for  obtaining  refined  range  infonnation  from  radar  by 
frequency  modulating  a  continuous  signal.  The  FMCW  modulation  is  readily  compatible 
with  a  wide  variety  of  solid-state  transmitters.  The  frequency  measurement  of  FMCW, 
which  is  perfonned  to  obtain  range  measurement,  can  be  performed  digitally  using  the 
FFT.  The  FMCW  signals  are  very  difficult  to  detect  with  conventional  intercept 
receivers.  [Ref.  14] 

The  frequency  modulation  used  by  the  radar  can  take  many  forms.  Linear  and 
sinusoidal  modulations  have  been  used  before.  Linear  frequency  modulation  is  very 
suitable  with  the  FFT  processor  to  get  the  range  information  of  targets.  The  linear  FMCW 
emitter  uses  a  continuous  100%  duty  cycle  waveform  so  that  both  the  target  range  and 
Doppler  information  can  be  measured.  The  FMCW  wavefonn  represents  the  best  use  of 
output  power  available  from  solid-state  devices.  Linear  FMCW  is  also  easier  to 
implement  than  phase  code  modulation  as  long  as  there  is  no  strict  demand  on  linearity 
specifications  over  the  modulation  bandwidth.  [Ref.  2] 

The  most  basic  Linear  FMCW  technique  varies  frequency  in  a  sawtooth  fashion, 
whereby  frequency  values  ramp  up  over  the  course  of  a  modulation  period  and  then  drop 
back  down  to  the  starting  frequency  to  begin  the  ramp  up  for  the  next  modulation  period. 
Several  repetitions  of  this  sequence  produce  a  sawtooth  shape  in  the  frequency  vs.  time 
plot.  By  processing  the  reflected  waveform,  one  can  extract  range  infonnation  about  the 
target.  [Ref.  15] 
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Figure  21  shows  a  functional  block  diagram  of  a  Triangular  FMCW  transmitter. 
A  sinusoidal  carrier  is  frequency-modulated  with  a  message  signal  that  takes  the 
triangular  shape  as  shown  in  Figure  22.  The  modulated  signal  is  then  corrupted  with 
channel  distortion,  which  is  assumed  to  be  additive  white  Gaussian  noise  (AWGN)  with 
zero  mean  and  variance. 


Sinusoidal 

Carrier 


Noise 


(AWGN) 


♦ 


FMCW 

Modulator 


*©■ 


>  Transmitted 
Waveform 


Triangle 

Generator 


Figure  2 1 .  Triangular  FMCW  Block  Diagram. 


Figure  22.  FMCW  Modulation  Signal  in  Time  Domain. 
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The  variables  in  the  diagram  are  as  follows. 

=  Carrier  frequency 
A/7  =  Modulation  bandwidth 
tm  =  Modulation  period 

The  up  ramp  and  down  ramp  of  each  triangle  must  be  constructed  separately  by  the 
triangle  generator.  The  modulating  waveform  for  the  up  ramp  can  be  modeled  as 


/«(*)=■ 


AF  AF 

fr - 1 - 1,  0 <t<t, 

2  t 

0,  elsewhere 


m 


(4.1) 


Similarly,  the  modulating  waveform  for  the  down  ramp  is  modeled  as 


fd(t)= 


,  AF  AF  _ 

/c+  ?  t,  0<t< tm 

^  hn 

0,  elsewhere 


(4.2) 


In  FM  systems,  the  transmit  phase  (instantaneous  frequency)  of  the  carrier  is  computed  as 


(f){t)  =  2n {  / (x)dx 
0 


(4.3) 


Now,  assuming  the  initial  conditions  of  (f>0  =  0  at  t  =  0  and  substituting  fu  (t)  and  frf(t) 
into  (4.3)  yields  the  RF  phases  for  the  up  ramp  and  down  ramp  as 


AF  AF  i 

<t>iXi)  =  2n\{fc-—  )-t+—r H,  0  <t<tn 


2 

AF 


AF 


4>d(t)  =  2nUfc+—  )-t~  —  t  C  0  <  t  <  t, 


(4.4) 


(4.5) 


It  is  then  a  simple  substitution  to  arrive  at  the  form  of  the  transmitted  signals  for 
the  up  and  down  ramps  [Ref.  14,  15] 
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(4.6) 


J„i(0  =  ^cos2^|(/c-^)-r  +  ^f2|, 
^i(0  =  ^cos2^-|(/c  +  ^)-r-^r|, 
^2  (0  =  ^4  sin  2tt  |(y;  -  r2 1 , 

{A  j7  a  f1  1 

(/c+— )-^-  — k 


0  <  t  <  tm 
0  <  t  <  tm 
0  <  t  <  tm 

o  <  t  <  tm 


(4.7) 

(4.8) 

(4.9) 


Like  those  codes  in  BPSK,  the  amplitude  A  is  equal  to  one.  The  most  important 
three  parameters  in  FMCW  code  generator  are  fc  ,  A F  and  tm  . 

Twelve  FMCW  signals  will  be  introduced  in  this  chapter.  They  are  numbered 
from  13  to  24  and  listed  in  Table  4.  The  “Mod.  BW”  means  the  modulation 
bandwidth;  250Hz  and  500Hz  are  used  in  this  table.  The  “Mod.  Period”  means  the 
modulation  period,  20  ms  and  30  ms  are  used  here.  A  1kHz  carrier  frequency  and 
7kHz  sampling  frequency  are  used  for  these  signals 


No. 

FMCW 

Mod.  BW  (Hz) 

Mod.  Period  (ms) 

SNR 

13 

F  1  7  250  20  s 

250 

20 

Signal  Only 

14 

F  1  7  250  20  0 

250 

20 

OdB 

15 

F  1  7  250  20  -6 

250 

20 

-6  dB 

16 

F  1  7  250  30  s 

250 

30 

Signal  Only 

17 

F  1  7  250  30  0 

250 

30 

OdB 

18 

F  1  7  250  30  -6 

250 

30 

-6  dB 

19 

F  1  7  500  20  s 

500 

20 

Signal  Only 

20 

F  1  7  500  20  0 

500 

20 

OdB 

21 

F  1  7  500  20  -6 

500 

20 

-6  dB 

22 

F  1  7  500  30  s 

500 

30 

Signal  Only 

23 

F  1  7  500  30  0 

500 

30 

OdB 

24 

F  1  7  500  30  -6 

500 

30 

-6  dB 

Table  4.  FMCW  Signals. 
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B.  WIGNER  DISTRIBUTION  FOR  FMCW 

Consider  the  same  twelve  FMCW  signals  in  Table  4  and  use  them  as  the  inputs 
for  the  Wigner  distribution  (WD).  The  WD  results  of  these  twelve  signals  are  shown  in 
Figure  23  to  26.  The  mesh  plots  show  the  frequency  domain  of  the  FMCW  signals  after 
the  WD.  The  contour  plots  show  both  the  frequency  domain  and  time  domain  of  the 
results. 

Take  Figure  23  as  an  example.  First,  the  carrier  frequency  fc  can  be  clearly 

found  by  the  location  of  the  highest  or  lowest  peak  value  in  Figure  23(a).  In  Figure 
23(b),  by  the  location  of  the  arc  shape  one  can  also  identify  the  carrier  frequency. 

Secondly,  the  modulation  bandwidth  A F  and  modulation  period  tm  have  scaled 
and  marked  in  Figures  23(b)  to  26(b) 
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Figure  23.  WD  for  FMCW  with  A F  =  250 Hz,tm  =  20 ms  Signal  Only  (a)  2D  Mesh 

Frequency  Domain,  (b)  Contour. 
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Figure  24.  WD  for  FMCW  with  A F  =  250 Hz,tm  =  30 ms  Signal  Only  (a)  2D  Mesh 

Frequency  Domain,  (b)  Contour. 
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Figure  25.  WD  for  FMCW  with  A F  =  500 Hz,tm  =  20 ms  Signal  Only  (a)  2D  Mesh 

Frequency  Domain,  (b)  Contour. 
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Figure  26.  WD  for  FMCW  with  A F  =  500 Hz,tm  =  30 ms  Signal  Only  (a)  2D  Mesh 

Frequency  Domain,  (b)  Contour. 
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C.  SUMMARY 

For  all  the  LPI  signals  in  this  thesis,  Wigner  Distribution  provides  better 
parameter  estimation  for  FMCW  than  the  others.  Without  any  other  analysis  beyond  the 
contour  plot,  one  can  see  the  most  two  important  parameters,  modulation  bandwidth  and 
modulation  period. 

Chapter  I  introduced  the  fact  that  the  Wigner  distribution  is  a  Time -Frequency 
Analysis  (TFA)  technique.  Since  the  FMCW  signals  have  the  character  that  both  the 
frequency  and  time  were  modulated,  using  the  Wigner  distribution  can  yield  excellent 
results  for  FMCW  signals,  even  at  SNR  equal  to  -6  dB.  See  the  technical  report  to  be 
published  soon  for  more  details.  [Ref.  16] 

Table  5  below  shows  the  detection  effectiveness  of  the  carrier  frequency, 
modulation  bandwidth  and  modulation  period  by  visible  measurement  from  Figure  23  to 
26.  They  are  done  by  the  visible  inspection  of  the  WD  processing  results. 

After  the  BPSK  and  FMCW  signals,  the  polyphase  signals  and  their  WD  analysis 
will  be  examined  in  the  next  chapter. 
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FMCW 

No. 

Carrier  Frequency 

Mod.  BW  (Hz) 

Mod.  Period  (ms) 

13 

100% 

100% 

100% 

14 

100% 

100% 

100% 

15 

100% 

90% 

100% 

16 

100% 

100% 

100% 

17 

100% 

100% 

100% 

18 

100% 

90% 

100% 

19 

100% 

100% 

100% 

20 

100% 

100% 

100% 

21 

100% 

90% 

100% 

22 

100% 

100% 

100% 

23 

100% 

100% 

100% 

24 

100% 

90% 

100% 

FMCW  effectiveness  on  detection 


150% 

130% 

110% 

<D 

O) 

CO 

90% 

c 

0) 

o 

70% 

a> 

CL 

50% 

30% 

10% 

□  Only  signal 
■  0  dB 

□  (-)  6  dB 


Carrier  freq.(Hz) 


Mod.Bandw  idth  (Hz) 

Parameters 


Mod.Period(ms) 


Table  5.  WD  Detection  Effectiveness  for  FMCW  Signals. 


51 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


52 


V.  POLYPHASE 


A.  POLYPHASE  CODE  MODULATION 

Polyphase  codes  have  many  useful  features  such  as  low  range-time  sidelobes, 
ease  of  implementation,  compatibility  with  digital  implementation,  low  cross-correlation 
between  codes,  the  Doppler  tolerance  of  the  frequency  modulation  codes  and 
compatibility  with  bandpass  limited  receivers.  The  largest  pulse  compression  ratio  that 
can  be  obtained  with  the  linear  Barker  codes  is  only  13  but  the  Polyphase  codes  are  not 
limited  to  any  finite  pulse  compression  ratio.  [Ref.  17] 

The  phase  of  the  subpulses  in  Polyphase  coded  pulse  compression  need  not  be 
restricted  to  the  two  levels  of  0  and  n  in  the  binary  phase  codes.  Polyphase  codes  produce 
lower  sidelobe  levels  than  the  binary  phase  codes  and  are  tolerant  to  Doppler  frequency 
shift  if  the  Doppler  frequencies  are  not  too  large.  [Ref.  4] 

The  Polyphase  code  waveforms  provide  a  class  of  frequency  derived  phase  coded 
waveforms  that  can  be  sampled  upon  reception  and  processed  digitally.  [Ref.  18] 
Polyphase  codes  also  addresses  the  reality  of  the  espionage  threat  as  well  as  the  electronic 
attack  threat.  Programmability  and  variety  are  two  code  characteristics  that  make 
espionage  a  more  difficult  task;  thus,  Polyphase  code  will  no  doubt  prove  to  be  the  LPI 
waveform  of  the  following  decades.  [Ref.  19] 

In  this  thesis,  the  polyphase  codes  including  Frank  code,  PI  code,  P2  code,  P3 
code,  and  P4  code  are  examined. 
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B. 


FRANK  CODE 


Frank  codes  are  a  family  of  polyphase  codes  that  are  closely  related  to  the  chirp 
and  have  been  used  successfully  in  LPI  radars.  The  Frank  code  has  a  length  N2 ,  which 
is  also  the  corresponding  pulse  compression  ratio.  If  i  is  the  number  of  the  sample  in  a 
given  frequency  and  j  is  the  number  of  the  frequency,  the  phase  of  the  z'th  sample  of  /th 
frequency  is  [Ref.  20] 

^=^(<■-1X7-1)  (5.1) 

where  i  =  1 :  N  and  j  =  1 :  N 


Twelve  Frank  code  signals  are  introduced  in  this  chapter.  They  are  numbered 
from  25  to  36  and  listed  in  Table  6.  In  this  chapter,  the  N  =  4  and  8  and  cycles  per  phase 
( CPP )  =  1  and  5  are  used.  A  1kHz  carrier  frequency  and  7kHz  sampling  frequency  are 
used  for  these  signals. 


No. 

FRANK 

N  (Phase) 

Cycles/Phase 

SNR 

25 

FR  1  7  4  1  s 

4 

1 

Signal  Only 

26 

FR  1  7  4  1  0 

4 

1 

OdB 

27 

FR  1  7  4  1  -6 

4 

1 

-6  dB 

28 

FR  1  7  4  5  s 

4 

5 

Signal  Only 

29 

FR  1  7  4  5  0 

4 

5 

OdB 

30 

FR  1  7  4  5  -6 

4 

5 

-6  dB 

31 

FR  1  7  8  1  s 

8 

1 

Signal  Only 

32 

FR  1  7  8  1  0 

8 

1 

OdB 

33 

FR  1  7  8  1  -6 

8 

1 

-6  dB 

34 

FR  1  7  8  5  s 

8 

5 

Signal  Only 

35 

FR  1  7  8  5  0 

8 

5 

OdB 

36 

FR  1  7  8  5  -6 

8 

5 

-6  dB 

Table  6.  Frank  Code  Signals. 
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1.  Wigner  Distribution  for  Frank  Code 

Consider  the  same  twelve  Frank  code  signals  in  Table  6  and  use  them  as  the 
inputs  for  Wigner  Distribution  (WD).  The  WD  results  of  these  twelve  signals  are  shown 
in  Figure  27  to  30.  The  mesh  plots  show  the  frequency  domain  of  the  Frank  code  signals 
after  WD.  The  contour  plots  show  both  the  frequency  domain  and  time  domain  of  the 
results. 


Take  Figure  27  as  an  example.  First,  the  carrier  frequency  fc  can  be  clearly 

found  by  the  location  of  the  highest  or  lowest  peak  value  in  Figure  27(a).  Secondly,  in 
Figure  27(b),  the  3dB  bandwidth  can  be  measured  as  1000Hz.  Since 


B  =  —  =  ^~ 
th  CPP 


(5.2) 


so,  the  cycles  per  phase  (CPP)  can  be  computed  by  CPP  = 


A 

B 


1000 

1000 


For  Frank  code 


9  N2 

T  =  N%=—  (5.3) 

since  the  code  period  ( T)  is  measured  as  0.016  in  Figure  27(b),  so  the  phase  length  (N2) 
can  be  calculated  from  N2  =  B-T  =  1000  •  0.016  =  16  .  The  slope  which  measured  in 
Figure  27(b)  as  1.6  x  10~6  sec2  can  be  also  use  to  determine  the  N2  and  CPP  by  the 
following  equation. 


T  _  N2  _  N2  CPP 2 

~b~1A~  J2 


(5.4) 
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Figure  27.  WD  for  Frank  Code  with  Phase  Length  =16,  CPP  =  1,  Signal  Only,  (No. 25) 
(a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Frequency  (Hz) 
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Figure  28.  WD  for  Frank  Code  with  Phase  Length  =  16,  CPP  =  5,  Signal  Only,  (No. 28) 
(a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  29.  WD  for  Frank  Code  with  Phase  Length  =  64,  CPP  =  1,  Signal  Only,  (No. 3 1) 
(a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  30.  WD  for  Frank  Code  with  Phase  Length  =  64,  CPP  =  5,  Signal  Only,  (No. 34) 
(a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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2.  Summary 


Applying  the  Wigner  Distribution,  one  can  quickly  know  the  carrier  frequency 
(/c)  of  the  Frank  codes.  The  carrier  frequency  is  always  the  first  parameter  need  to  be 
identified.  Once  the  carrier  frequency  is  found,  by  the  bandwidth  can  detennine  the  CPP 
;  by  the  code  period  or  the  slope  can  determine  the  phase  length(  N2 ). 


In  Figure  28(b)  (A7  =16,  CPP= 5)  only  bandwidth  is  obtained  but  not  the  code 
period.  The  reason  is  that  the  sample  number  ( #sample )  is  not  enough.  From  (5.3),  the 

code  period  is  T  =  ^—  =  -'-^-  =  0.080  sec  ,  hence  only  when  ^samP^e  >  j  the  code 
B  200  fs 


#samvle  512 

period  can  be  shown  in  the  figure.  In  this  chapter,  the  — — — —  =  — —  =  0.073  <  0.080 . 


/, 


7000 


Table  7  below  shows  the  detection  effectiveness  of  the  carrier  frequency,  3dB 
bandwidth,  code  period  and  phase  length  from  Figure  27  to  30.  They  are  done  by  the 
visible  inspection  of  the  WD  results. 
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Frank 

No. 

Carrier  Frequency 

Bandwidth 

Code  Period 

Phase  Fength 

25 

100% 

100% 

100% 

100% 

26 

90% 

110% 

100% 

100% 

27 

110% 

110% 

100% 

100% 

28 

95% 

100% 

0% 

0% 

29 

98% 

95% 

0% 

0% 

30 

0% 

0% 

0% 

0% 

31 

100% 

100% 

100% 

100% 

32 

95% 

105% 

100% 

100% 

33 

98% 

105% 

100% 

100% 

34 

100% 

100% 

0% 

0% 

35 

100% 

0% 

0% 

0% 

36 

100% 

0% 

0% 

0% 

Frank  code  effectiveness  on  detection 


Parameters 


□  Only  signal 
■  0  dB 

□  (-)  6  dB 


Table  7.  WD  Detection  Effectiveness  for  Frank  Code  Signals. 
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c. 


PI  CODE 


In  case  of  a  double  sideband  detection  (local  oscillator  is  at  band  center)  of  a  step 
approximation  of  a  linear  frequency  modulation,  which  is  the  PI  code.  The  PI  code  also 
consists  of  N2  elements.  If  i  is  the  number  of  the  sample  in  a  given  frequency  and  j  is 
the  number  of  the  frequency,  the  phase  of  the  /th  sample  of y'th  frequency  is  [Ref.  20] 

Aj  =^-[N-(2j-mU-l)N  +  (i~l)\  (5-5) 

where  i  =  1 :  N  and  j  =  1 :  N 


Twelve  PI  code  signals  are  introduced  in  this  chapter.  They  are  numbered  from 
37  to  48  and  listed  in  Table  8.  A  1kHz  carrier  frequency  and  7kHz  sampling  frequency 
are  used  for  these  signals. 


No. 

PI 

N  (Phase) 

Cycles/Phase 

SNR 

37 

PI  1  7  4  1  s 

4 

1 

Signal  Only 

38 

PI  1  7  4  1  0 

4 

1 

OdB 

39 

PI  1  7  4  1  -6 

4 

1 

-6  dB 

40 

PI  1  7  4  5  s 

4 

5 

Signal  Only 

41 

PI  1  7  4  5  0 

4 

5 

OdB 

42 

PI  1  7  4  5  -6 

4 

5 

-6  dB 

43 

PI  1  7  8  1  s 

8 

1 

Signal  Only 

44 

PI  1  7  8  1  0 

8 

1 

OdB 

45 

PI  1  7  8  1  -6 

8 

1 

-6  dB 

46 

PI  1  7  8  5  s 

8 

5 

Signal  Only 

47 

PI  1  7  8  5  0 

8 

5 

OdB 

48 

PI  1  7  8  5  -6 

8 

5 

-6  dB 

Table  8.  PI  Code  Signals. 
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1.  Wigner  Distribution  for  PI  Code 


Consider  the  same  twelve  PI  code  signals  in  Table  8  and  use  them  as  the  inputs 
for  Wigner  Distribution  (WD).  The  WD  results  of  these  twelve  signals  are  showed  in 
Figure  31  to  34.  The  mesh  plots  show  the  frequency  domain  of  the  PI  code  signals  after 
WD.  The  contour  plots  show  both  the  frequency  domain  and  time  domain  of  the  results. 


Take  Figure  31  as  an  example.  First,  the  carrier  frequency  fc  can  be  clearly 

found  by  the  location  of  the  highest  or  lowest  peak  value  in  Figure  31(a).  Secondly,  in 
Figure  3 1(b),  the  3dB  bandwidth  can  be  measured  as  1000Hz.  Since 


B  =  —  =  ^~ 
th  CPP 


(5.6) 


so,  the  cycles  per  phase  (CPP)  can  be  computed  by  CPP  = 


A 

B 


1000 

1000 


For  PI  code 


,  N 2 

T  =  N%=—  (5.7) 

since  the  code  period  ( T)  is  measured  as  0.016  in  Figure  31(b),  so  the  phase  length  (N2) 
can  be  calculated  from  N2  =  B-T  =  1000 •  0.016  =  16  .  The  slope  which  measured  in 
Figure  31(b)  as  1.6xl0~6sec2  can  be  also  use  to  determine  the  N2  and  CPP  by  the 
following  equation. 

C_T  N 2  N2  -  CPP2 

S~1  =  W~  f;  <58) 
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Figure  31.  WD  for  PI  Code  with  Phase  Length  =16,  CPP  =  1,  Signal  Only,  (No. 37)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  32.  WD  for  PI  Code  with  Phase  Length  =  16,  CPP  =  5,  Signal  Only,  (No.40)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  33.  WD  for  PI  Code  with  Phase  Length  =  64,  CPP  =  1,  Signal  Only,  (No.43)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  34.  WD  for  PI  Code  with  Phase  Length  =  64,  CPP  =  5,  Signal  Only,  (No.46)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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2.  Summary 


Applying  the  Wigner  Distribution,  one  can  quickly  know  the  carrier  frequency 
( fc )  of  the  Pi  codes.  The  carrier  frequency  is  always  the  first  parameter  need  to  be 
identified.  Once  the  carrier  frequency  is  found,  by  the  bandwidth  can  detennine  the  CPP 
;  by  the  code  period  or  the  slope  can  determine  the  phase  length(  N2 ). 


In  Figure  32(b)  (A7  =16,  CPP= 5)  only  bandwidth  is  obtained  but  not  the  code 
period.  The  reason  is  that  the  sample  number  ( #sample )  is  not  enough.  From  (5.7),  the 

code  period  is  T  =  ^—  =  -'-^-  =  0.080  sec  ,  hence  only  when  ^samP^e  >  j  the  code 
B  200  fs 


#samvle  512 

period  can  be  shown  in  the  figure.  In  this  chapter,  the  — — — —  =  — —  =  0.073  <  0.080 . 


/, 


7000 


Table  9  below  shows  the  detection  effectiveness  of  the  carrier  frequency,  3dB 
bandwidth,  code  period  and  phase  length  from  Figure  31  to  34.  They  are  done  by  the 
visible  inspection  of  the  WD  results. 
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PI 

No. 

Carrier  Frequency 

Bandwidth 

Code  Period 

Phase  Length 

37 

120% 

70% 

100% 

100% 

38 

110% 

80% 

100% 

100% 

39 

0% 

0% 

0% 

0% 

40 

100% 

100% 

0% 

0% 

41 

100% 

90% 

0% 

0% 

42 

0% 

0% 

0% 

0% 

43 

90% 

50% 

100% 

100% 

44 

90% 

60% 

100% 

100% 

45 

0% 

0% 

0% 

0% 

46 

100% 

100% 

0% 

0% 

47 

100% 

120% 

0% 

0% 

48 

0% 

0% 

0% 

0% 

PI  code  effectiveness  on  detection 


Parameters 


□  Only  signal 

□  0  dB 

□  (-)  6  dB 


Table  9.  WD  Detection  Effectiveness  for  PI  Code  Signals. 
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D.  P2  CODE 


The  P2  code  also  consists  N2  elements.  If  i  is  the  number  of  the  sample  in  a 
given  frequency  and  j  is  the  number  of  the  frequency,  the  phase  of  the  /th  sample  of  /th 
frequency  is  [Ref.  20] 


(5.9) 


where  i  =  1 :  N  and  j  =  1 :  N 


Twelve  P2  code  signals  are  introduced  in  this  chapter.  They  are  numbered  from 


49  to  60  and  listed  in  Table  10.  Phase  length  N2=  16,  64  and  CPP=  1,  5  are  used  in  this 
chapter.  A  1kHz  carrier  frequency  and  7kHz  sampling  frequency  are  used  for  these 
signals. 


No. 

P2 

N  (Phase) 

Cycles/Phase 

SNR 

49 

P2  1  7  4  1  s 

4 

1 

Signal  Only 

50 

P2  1  7  4  1  0 

4 

1 

OdB 

51 

P2  1  7  4  1  -6 

4 

1 

-6  dB 

52 

P2  1  7  4  5  s 

4 

5 

Signal  Only 

53 

P2  1  7  4  5  0 

4 

5 

OdB 

54 

P2  1  7  4  5  -6 

4 

5 

-6  dB 

55 

P2  1  7  8  1  s 

8 

1 

Signal  Only 

56 

P2  1  7  8  1  0 

8 

1 

OdB 

57 

P2  1  7  8  1  -6 

8 

1 

-6  dB 

58 

P2  1  7  8  5  s 

8 

5 

Signal  Only 

59 

P2  1  7  8  5  0 

8 

5 

OdB 

60 

P2  1  7  8  5  -6 

8 

5 

-6  dB 

Table  10.  P2  Code  Signals. 
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1.  Wigner  Distribution  for  P2  Code 


Consider  the  same  twelve  P2  code  signals  in  Table  10  and  use  them  as  the  inputs 
for  Wigner  Distribution  (WD).  The  WD  results  of  these  twelve  signals  are  showed  in 
Figure  35  to  38.  The  mesh  plots  show  the  frequency  domain  of  the  P2  code  signals  after 
WD.  The  contour  plots  show  both  the  frequency  domain  and  time  domain  of  the  results. 

Take  Figure  35  as  an  example.  First,  the  carrier  frequency  fc  can  be  clearly 

found  by  the  location  of  the  highest  or  lowest  peak  value  in  Figure  35(a).  Secondly,  in 
Figure  35(b),  the  3dB  bandwidth  can  be  measured  as  1000Hz.  Since 

1  f 

B  =  -  =  ^~  (5.10) 

tb  CPP 

so,  the  cycles  per  phase  (CPP)  can  be  computed  by  CPP  =  =  ~~  =  1  For  P2  code 


,  N 2 

T  =  N\=—  (5.11) 

since  the  code  period  ( T)  is  measured  as  0.016  in  Figure  35(b),  so  the  phase  length  (N2) 
can  be  calculated  from  N2  =  B-T  =  1000 •  0.016  =  16  .  The  slope  which  measured  in 
Figure  35(b)  as  1.6x  10~6  sec2  can  be  also  use  to  determine  the  N2  and  CPP  by  the 
following  equation. 


T  N2  _N2- CPP 2 


(5.12) 
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iP) 

Figure  35.  WD  for  P2  Code  with  Phase  Length  =16,  CPP  =  1,  Signal  Only,  (No.49)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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meshofP21745s 
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contour  of  P21745s 
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Figure  36.  WD  for  P2  Code  with  Phase  Length  =  16,  CPP  =  5,  Signal  Only,  (No. 52)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 


73 


500 


1000 


007 


1500  2000 

Frequency  (Hz) 

(a) 

contour  of  P2, 


2500 


3000 


1781s 


Bandwidth  ( B )  =  1 000  ]  iz 


(T)  =-t).064  S6c 


f  .Tt  ** 

Slop</(S)  =  -*  =  64  x  F0_fi  sec 
B 


800 


600 


§>200 

03 


-200 


meshofP21781s 


0  500  1000  1500  2000  2500  3000 

Frequency  (Hz) 

(b) 

Figure  37.  WD  for  P2  Code  with  Phase  Length  =  64,  CPP  =  1,  Signal  Only,  (No. 55)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  38.  WD  for  P2  Code  with  Phase  Length  =  64,  CPP  =  5,  Signal  Only,  (No. 58)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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2.  Summary 


Applying  the  Wigner  Distribution,  one  can  quickly  know  the  carrier  frequency 
( fc )  of  the  P2  codes.  The  carrier  frequency  is  always  the  first  parameter  need  to  be 
identified.  Once  the  carrier  frequency  is  found,  by  the  bandwidth  can  detennine  the  CPP 
;  by  the  code  period  or  the  slope  can  determine  the  phase  length(  N2 ). 


In  Figure  36(b)  (A7  =16,  CPP= 5)  only  bandwidth  is  obtained  but  not  the  code 
period.  The  reason  is  that  the  sample  number  ( #sample )  is  not  enough.  From  (5.11),  the 

code  period  is  T  =  ^—  =  -'-^-  =  0.080  sec  ,  hence  only  when  ^samP^e  >  j  the  code 
B  200  fs 


#samvle  512 

period  can  be  shown  in  the  figure.  In  this  chapter,  the  — — — —  =  — —  =  0.073  <  0.080 . 


7000 


Table  1 1  below  shows  the  detection  effectiveness  of  the  carrier  frequency,  3dB 
bandwidth,  code  period  and  phase  length  from  Figure  35  to  38.  They  are  done  by  the 
visible  inspection  of  the  WD  results. 
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P2 

No. 

Carrier  Frequency 

Bandwidth 

Code  Period 

Phase  Length 

49 

100% 

70% 

100% 

100% 

50 

110% 

60% 

90% 

100% 

51 

0% 

0% 

0% 

0% 

52 

100% 

100% 

0% 

0% 

53 

100% 

0% 

0% 

0% 

54 

95% 

0% 

0% 

0% 

55 

100% 

100% 

100% 

100% 

56 

100% 

110% 

110% 

100% 

57 

0% 

110% 

110% 

100% 

58 

100% 

100% 

0% 

0% 

59 

100% 

90% 

0% 

0% 

60 

120% 

0% 

0% 

0% 

P2  code  effectiveness  on  detection 


Parameters 


□  Only  signal 

□  0  dB 

□  (-)  6  dB 


Table  11.  WD  Detection  Effectiveness  for  P2  Code  Signals. 
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E. 


P3  CODE 


In  the  case  of  linear  frequency  modulation  waveforms,  the  conceptual  coherent 
detection  using  single  sideband  detection  and  sampling  at  the  Nyquist  rate  yield  a 
polyphase  code  call  the  P3.  The  phase  of  the  zth  sample  of  the  P3  code  is  given  by  [Ref. 
20] 


4=^0'— l)(i-l) 

where  i  =  1 :  N 


(5.13) 


Twelve  P3  code  signals  are  introduced  in  this  chapter.  They  are  numbered  from 
61  to  72  and  listed  in  Table  12.  In  this  chapter,  the  phase  length  N  =  16  and  64  and 
cycles  per  phase  ( CPP )  =  1  and  5  are  used.  A  1kHz  carrier  frequency  and  7kHz  sampling 
frequency  are  used  for  these  signals. 


No. 

P3 

N  (Phase) 

Cycles/Phase 

SNR 

61 

P3  1  7  16  1  s 

16 

1 

Signal  Only 

62 

P3  1  7  16  1  0 

16 

1 

OdB 

63 

P3  1  7  16  1  -6 

16 

1 

-6  dB 

64 

P3  1  7  16  5  s 

16 

5 

Signal  Only 

65 

P3  1  7  16  5  0 

16 

5 

OdB 

66 

P3  1  7  16  5  -6 

16 

5 

-6  dB 

67 

P3  1  7  64  1  s 

64 

1 

Signal  Only 

68 

P3  1  7  64  1  0 

64 

1 

OdB 

69 

P3  1  7  64  1  -6 

64 

1 

-6  dB 

70 

P3  1  7  64  5  s 

64 

5 

Signal  Only 

71 

P3  1  7  64  5  0 

64 

5 

OdB 

72 

P3  1  7  64  5  -6 

64 

5 

-6  dB 

Table  12.  P3  Code  Signals. 
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1.  Wigner  Distribution  for  P3  Code 


Consider  the  same  twelve  P3  code  signals  in  Table  12  and  use  them  as  the  inputs 
for  Wigner  Distribution  (WD).  The  WD  results  of  these  twelve  signals  are  showed  in 
Figure  39  to  42.  The  mesh  plots  show  the  frequency  domain  of  the  P3  code  signals  after 
WD.  The  contour  plots  show  both  the  frequency  domain  and  time  domain  of  the  results. 

Take  Figure  39  as  an  example.  First,  the  carrier  frequency  fc  can  be  clearly 
found  by  the  location  of  the  highest  or  lowest  peak  value  in  Figure  39(a).  Secondly,  in 
Figure  39(b),  the  3dB  bandwidth  can  be  measured  as  1000Hz.  Since 

1  f 

B  =  —  =  ^~  (5.14) 

tb  CPP 

so,  the  cycles  per  phase  (CPP)  can  be  computed  by  CPP  =  =  ~~  =  1  For  P3  code 


T  =  Nt„=j  (5-15) 

since  the  code  period  (7)  is  measured  as  0.016  in  Figure  39(b),  so  the  phase  length  (N) 
can  be  calculated  from  N  =  B  -T  =  1000 •  0.016  =  16  .  The  slope  which  measured  in 
Figure  39(b)  as  1.6  x  10~6  sec2  can  be  also  use  to  determine  the  N2  and  CPP  by  the 
following  equation. 


T  _  N  _  N  ■  CPP  2 
B~BY~~ff~ 


(5.16) 
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Figure  39.  WD  for  P3  Code  with  Phase  Length  =  16,  CPP  =  1,  Signal  Only,  (No. 61)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  40.  WD  for  P3  Code  with  Phase  Length  =  16,  CPP  =  5,  Signal  Only,  (No. 64)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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mesh  of  P3.,  7g41 
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Figure  41.  WD  for  P3  Code  with  Phase  Length  =  64,  CPP  =  1,  Signal  Only,  (No. 67)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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mesh  of  P317g45s 


(a) 

contour  of  P317645 


(b) 

Figure  42.  WD  for  P3  Code  with  Phase  Length  =  64,  CPP  =  5,  Signal  Only,  (No. 70)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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2.  Summary 


Applying  the  Wigner  Distribution,  one  can  quickly  know  the  carrier  frequency 
( fc )  of  the  P3  codes.  The  carrier  frequency  is  always  the  first  parameter  need  to  be 

identified.  Once  the  carrier  frequency  is  found,  by  the  bandwidth  can  detennine  the  CPP 
;  by  the  code  period  or  the  slope  can  determine  the  phase  length(  N  ). 


In  Figure  40(b)  (A  =16,  CPP=5)  only  bandwidth  is  obtained  but  not  the  code 
period.  The  reason  is  that  the  sample  number  ( #sample )  is  not  enough.  From  (5.15),  the 

code  period  is  T  =  —  =  =  0.080  sec ,  hence  only  when  —  >  T  ,  the  code  period 

B  200  fs 


O'  scuTivlc  512 

can  be  shown  in  the  figure.  In  this  chapter,  the  — - = - =  0.073  <  0.080  . 


fs 


7000 


Table  13  below  shows  the  detection  effectiveness  of  the  carrier  frequency,  3dB 
bandwidth,  code  period  and  phase  length  from  Figure  39  to  42.  They  are  done  by  the 
visible  inspection  of  the  WD  results. 
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P3 

No. 

Carrier  Frequency 

Bandwidth 

Code  Period 

Phase  Length 

61 

95% 

100% 

100% 

100% 

62 

95% 

110% 

100% 

100% 

63 

80% 

110% 

110% 

100% 

64 

100% 

100% 

0% 

0% 

65 

100% 

120% 

0% 

0% 

66 

90% 

0% 

0% 

0% 

67 

110% 

100% 

90% 

100% 

68 

110% 

100% 

90% 

100% 

69 

120% 

110% 

90% 

100% 

70 

100% 

100% 

0% 

0% 

71 

100% 

0% 

0% 

0% 

72 

100% 

0% 

0% 

0% 

P3  code  effectiveness  on  detection 


Parameters 


□  Only  signal 
■  0  dB 

□  (-)  6  dB 


Table  13.  WD  Detection  Effectiveness  for  P3  Code  Signals. 
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F.  P4  CODE 

The  P4  code  consists  of  the  discrete  phase  of  the  linear  chirp  waveform  taken  at 
specific  time  intervals  and  exhibits  the  same  range  Doppler  coupling  associated  with  the 
chirp  waveform.  The  peak  sidelobe  levels  are  lower  than  those  of  the  unweighted  chirp 
waveform.  The  code  element  of  the  P4  code  are  given  by  [Ref.  20] 

1)  (5.17) 

where  i  =  1 :  N 


Twelve  P4  code  signals  are  introduced  in  this  chapter.  They  are  numbered  from 
73  to  84  and  listed  in  Table  14.  The  phase  length  or  phase  elements  N  are  also  16  or  64; 
CPP  is  1  or  5  for  P4  codes  in  this  chapter.  A  1kHz  carrier  frequency  and  7kHz  sampling 
frequency  are  used  for  these  signals. 


No. 

P4 

N  (Phase) 

Cycles/Phase 

SNR 

73 

P4  1  7  16  1  s 

16 

1 

Signal  Only 

74 

P4  1  7  16  1  0 

16 

1 

OdB 

75 

P4  1  7  16  1  -6 

16 

1 

-6  dB 

76 

P4  1  7  16  5  s 

16 

5 

Signal  Only 

77 

P4  1  7  16  5  0 

16 

5 

OdB 

78 

P4  1  7  16  5  -6 

16 

5 

-6  dB 

79 

P4  1  7  64  1  s 

64 

1 

Signal  Only 

80 

P4  1  7  64  1  0 

64 

1 

OdB 

81 

P4  i  7  64  1  -6 

64 

1 

-6  dB 

82 

P4  1  7  64  5  s 

64 

5 

Signal  Only 

83 

P4  1  7  64  5  0 

64 

5 

OdB 

84 

P4  1  7  64  5  -6 

64 

5 

-6  dB 

Table  14.  P4  Code  Signals. 
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1.  Wigner  Distribution  for  P4  Code 


Consider  the  same  twelve  P4  code  signals  in  Table  14  and  use  them  as  the  inputs. 
The  WD  results  of  these  twelve  signals  are  showed  in  Figure  43  to  46.  The  mesh  plots 
show  the  frequency  domain  of  the  P4  code  signals  after  WD.  The  contour  plots  show 
both  the  frequency  domain  and  time  domain  of  the  results. 

Take  Figure  43  as  an  example.  First,  the  carrier  frequency  fc  can  be  clearly 

found  by  the  location  of  the  highest  or  lowest  peak  value  in  Figure  43(a).  Secondly,  in 
Figure  43  (b),  the  3dB  bandwidth  can  be  measured  as  1000Hz.  Since 

1  f 

B  =  -  =  ^~  (5.18) 

tb  CPP 

so,  the  cycles  per  phase  (CPP)  can  be  computed  by  CPP  =  =  =  '  F°r  ^  code 


T  =  Ntb  = 


N_ 

B 


(5.19) 


since  the  code  period  ( T)  is  measured  as  0.016  in  Figure  43  (b),  so  the  phase  length  (  N  ) 
can  be  calculated  from  N  =  B  -T  =  1000 •  0.016  =  16  .  The  slope  which  measured  in 
Figure  43  (b)  as  1.6  x  10~6  sec2  can  be  also  use  to  determine  the  N2  and  CPP  by  the 
following  equation. 


T  _  N  _  N  ■  CPP  2 
B~BY~~ff~ 


(5.20) 
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mesh  of  P417161s 


(a) 


contour  of  P417161$ 


0) 

Figure  43.  WD  for  P4  Code  with  Phase  Length  =16,  CPP  =  1,  Signal  Only,  (No. 73)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 


88 


mesh  of  P417165s 
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Figure  44.  WD  for  P4  Code  with  Phase  Length  =  16,  CPP  =  5,  Signal  Only,  (No. 76)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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mesh  ofP41/641s 


(a) 

contour  of  P417641s 


Frequency  (Hz) 


0) 

Figure  45.  WD  for  P4  Code  with  Phase  Length  =  64,  CPP  =  1,  Signal  Only,  (No. 79)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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mesh  ofP417645s 
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Figure  46.  WD  for  P4  Code  with  Phase  Length  =  64,  CPP  =  5,  Signal  Only,  (No. 82)  (a) 

2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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2.  Summary 


Applying  the  Wigner  Distribution,  one  can  quickly  know  the  carrier  frequency 
( fc )  of  the  P4  codes.  The  carrier  frequency  is  always  the  first  parameter  need  to  be 

identified.  Once  the  carrier  frequency  is  found,  by  the  bandwidth  can  detennine  the  CPP 
;  by  the  code  period  or  the  slope  can  determine  the  phase  length(  N  ). 


In  Figure  44(b)  (A  =16,  CPP=5)  only  bandwidth  is  obtained  but  not  the  code 
period.  The  reason  is  that  the  sample  number  ( #sample )  is  not  enough.  From  (5.19),  the 

code  period  is  T  =  —  =  =  0.080  sec ,  hence  only  when  —  >  T  ,  the  code  period 

B  200  fs 


O'  scuTivlc  512 

can  be  shown  in  the  figure.  In  this  chapter,  the  — - = - =  0.073  <  0.080  . 
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Table  9  below  shows  the  detection  effectiveness  of  the  carrier  frequency,  3dB 
bandwidth,  code  period  and  phase  length  from  Figure  43  to  46.  They  are  done  by  the 
visible  inspection  of  the  WD  results. 
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P4 

No. 

Carrier  Frequency 

Bandwidth 

Code  Period 

Phase  Length 

73 

90% 

100% 

100% 

100% 

74 

90% 

95% 

110% 

100% 

75 

80% 

90% 

110% 

100% 

76 

100% 

100% 

0% 

0% 

77 

95% 

0% 

0% 

0% 

78 

100% 

0% 

0% 

0% 

79 

110% 

100% 

100% 

100% 

80 

115% 

105% 

95% 

100% 

81 

125% 

110% 

90% 

100% 

82 

100% 

100% 

0% 

0% 

83 

100% 

110% 

0% 

0% 

84 

100% 

120% 

0% 

0% 

P4  code  effectiveness  on  detection 


□  Only  signal 
■  0  dB 

□  (-)  6  dB 


Parameters 


Table  15.  WD  Detection  Effectiveness  for  P4  Code  Signals. 
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G. 


COMPARISON  OF  POLYPHASE  CODES 


Table  16  shows  the  phases  for  all  the  different  polyphase  codes  discussed  in  this 

thesis. 


Frank  Code 

2  7T 

fa  j  =  (/  - !)(  /  - 1)  where  i  =  1 :  N  and  j  =  1 :  N 

PI  Code 

^ +  </-!)]  »’here  i-l-.N  end  J-l-.N 

P2  Code 

faij  =^^(2i-l-N){2j -\-N)  where  i  =  \:N  and  j  =  l:N 

P3  Code 

6  =  —  (i  —  l)(i  - 1)  where  i  =  1 :  N 

N 

P4  Code 

(  '  1 

a ).  =  n  - 7i  (i  - 1)  where  i  =  1 :  N 

N 

Table  16.  Phases  of  All  the  Polyphase  Code. 


Figure  47  compares  the  WD  of  all  the  polyphase  codes  with  same  CPP  =  1  and 
phase  length  =16.  Figure  48  is  the  zoom  in  figure  for  Figure  47.  Figure  49  compares  the 
WD  of  all  the  polyphase  codes  with  same  CPP  =  5  and  phase  length  =16.  Figure  50  is 
the  zoom  in  figure  for  Figure  49.  Figure  51  compares  the  WD  of  all  the  polyphase  codes 
with  same  CPP  =  1  and  phase  length  =  64.  Figure  52  is  the  zoom  in  figure  for  Figure  5 1 . 
Figure  53  compares  the  WD  of  all  the  polyphase  codes  with  same  CPP  =  5  and  phase 
length  =  64.  Figure  54  is  the  zoom  in  figure  for  Figure  53. 
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Figure  47.  Wigner  Distribution  for  Polyphase  Codes  with  CPP  =  1  and  Phase  Length  = 
16  (a)  Frank  Code,  (b)  PI  Code,  (c)  P2  Code,  (d)  P3  Code,  (e)  P4  Code. 
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Figure  48.  Zoom  In  for  Figure  47  (a)  Frank  Code,  (b)  PI  Code,  (c)  P2  Code,  (d)  P3  Code. 

(e)  P4  Code. 
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Figure  49.  Wigner  Distribution  for  Polyphase  Codes  with  CPP  =  5  and  Phase  Length  = 
16  (a)  Frank  Code,  (b)  PI  Code,  (c)  P2  Code,  (d)  P3  Code,  (e)  P4  Code. 
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Figure  50.  Zoom  In  for  Figure  49  (a)  Frank  Code,  (b)  PI  Code,  (c)  P2  Code,  (d)  P3  Code. 

(e)  P4  Code. 
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Figure  5 1 .  Wigner  Distribution  for  Polyphase  Codes  with  CPP  =  1  and  Phase  Length 
64  (a)  Frank  Code,  (b)  PI  Code,  (c)  P2  Code,  (d)  P3  Code,  (e)  P4  Code. 
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Figure  52.  Zoom  In  for  Figure  51  (a)  Frank  Code,  (b)  PI  Code,  (c)  P2  Code,  (d)  P3  Code. 

(e)  P4  Code. 
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Figure  53.  Wigner  Distribution  for  Polyphase  Codes  with  CPP  =  5  and  Phase  Length 
64  (a)  Frank  Code,  (b)  PI  Code,  (c)  P2  Code,  (d)  P3  Code,  (e)  P4  Code. 
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Figure  54.  Zoom  In  for  Figure  53  (a)  Frank  Code,  (b)  PI  Code,  (c)  P2  Code,  (d)  P3  Code. 

(e)  P4  Code. 
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Some  features  are  observed  from  the  above  eight  figures.  First,  for  all  the 
polyphase  signals  with  CPP  =1  (Figures  47,  48,  51  and  52),  the  P2  codes  have  negative 
slope,  which  is  different  from  other  polyphase  signals.  Secondly,  for  the  WD  analysis, 
when  CPP  =1  (Figures  47,  48,  51  and  52),  the  PI  code  is  the  most  difficult  to  determine 
the  3  dB  bandwidth  as  well  as  the  code  period.  On  the  contrary,  the  P3  and  P4  code  can 
provide  very  clear  parameters  in  Figure  48  and  52.  Thirdly,  when  comparing  with  all  the 
zoom  in  figures  (Figures  48,  50,  52  and  54),  the  P3  code  and  the  P4  code  have  the  best 
parameter  resolution  among  the  live  polyphase  signals. 
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VI.  COSTAS  FREQUENCY  HOPPING 


A.  COSTAS 

A  Costas  array  is  an  nxn  array  of  frequencies  and  times  with  exactly  one 
frequency  in  each  row  and  column,  and  with  distinct  vector  differences  between  all  pairs 
of  frequencies.  As  a  frequency-hopping  pattern  for  radar  or  sonar,  a  Costas  array  has  an 
optimum  ambiguity  function,  since  any  translation  of  the  array  parallel  to  the  coordinate 
axes  produces  at  most  one  out  of  phase  coincidence.  [Ref.  21] 

In  a  frequency  hopping  system,  the  signal  consists  of  one  or  more  frequencies 
being  chosen  from  a  set  {fl,f2,—fm}  of  available  frequencies,  for  transmission  at  each 

of  a  set  {t1,t2,....tn}  of  consecutive  time  intervals.  For  modeling  purposes,  it  is  reasonable 
to  consider  the  situation  in  which  m  =  n,  and  a  different  one  of  n  equally  spaced 
frequencies  is  transmitted  during  each  of  the  n  equal  duration  time  intervals 

{tx,t2,....tn}  .  Such  a  signal  is  represented  by  a  nxn  permutation  matrix  A,  where  the  n 
rows  correspond  to  the  n  frequencies,  the  n  columns  correspond  to  the  n  intervals,  and  the 
entry  atj  equals  1  means  transmission  and  0  otherwise. 

For  example  if  the  coding  sequence  {  aj  }  is  {4,7, 1,6, 5, 2, 3,},  then  the  coding 

matrix,  the  difference  matrix  and  the  ambiguity  sidelobe  matrix  of  a  Costas  signal  are 
shown  as  follows.  [Ref.  22] 
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Figure  55.  The  Coding  Matrix,  Difference  Matrix  and  Ambiguity  Sidelobe  Matrix  of  a 

Costas  Signal. 


The  element  of  the  difference  matrix  in  row  i  and  column  j,  is 


Di.i=ai-j-aP  i  +  j^N  (6.1) 

with  blanks  at  the  remaining  locations.  Thus  the  first  row  is  formed  by  taking  differences 
between  adjacent  terms  in  the  sequence,  the  second  row  by  taking  differences  between 
next  adjacent  terms,  and  so  on.  The  corresponding  values  in  the  difference  matrix  and 
sidelobe  matrix  were  marker  with  different  color  circles.  [Ref.  22] 


The  difference  matrix  not  only  detennines  whether  a  sequence  is  Costas  or  not, 
but  it  also  identifies  the  major  sidelobe  locations  in  the  positive  delay  slots  of  the 
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ambiguity  function.  The  sidelobes  in  the  negative  delay  slots  are  obtained  by  applying 


the  rule  of  symmetry  with  respect  with  the  origin.  [Ref.  22] 


Twelve  Costas  signals  are  introduced  in  this  chapter.  They  are  numbered  from  85 
to  96  and  listed  in  Table  17.  The  “Sequence”  means  the  frequency  hopping  sequence. 
There  are  two  different  sequences  used  in  this  table.  They  are  the  carrier  frequency  in 
order  of  the  number  in  the  sequence.  The  unit  of  these  carrier  frequencies  is  kHz.  The 
“Cycles/  Frequency”  means  the  number  of  cycles  per  frequency,  the  number  10  and  20 
are  used  here.  The  sampling  frequency  for  these  Costas  signals  is  15kHz  for  sequence 
4716523  and  17kHz  for  sequence  2638751. 


No. 

COSTAS 

Sequence 

Cy  c  les/F  requency 

SNR 

85 

C  1  15  10  s 

4716523 

10 

Signal  Only 

86 

C  1  15  10  0 

4716523 

10 

OdB 

87 

C  1  15  10  -6 

4716523 

10 

-6  dB 

88 

C  1  15  20  s 

4716523 

20 

Signal  Only 

89 

C  1  15  20  0 

4716523 

20 

OdB 

90 

C  1  15  20  -6 

4716523 

20 

-6  dB 

91 

C  2  17  10  s 

2638751 

10 

Signal  Only 

92 

C  2  17  10  0 

2638751 

10 

OdB 

93 

C  2  17  10  -6 

2638751 

10 

-6  dB 

94 

C  2  17  20  s 

2638751 

20 

Signal  Only 

95 

C  2  17  20  0 

2638751 

20 

OdB 

96 

C  2  17  20  -6 

2638751 

20 

-6  dB 

Table  17.  COSTAS  Signals. 
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B. 


WIGNER  DISTRIBUTION  FOR  COSTAS 


Consider  the  same  twelve  Costas  frequency-hopping  signals  in  Table  17  and  use 
them  as  the  inputs  for  Wigner  Distribution  (WD).  The  WD  results  of  signals  number  85, 
88,  91,  94  are  showed  in  Figures  56  to  59.  The  mesh  plots  show  the  frequency  domain  of 
the  Costas  frequency-hopping  signals  after  WD.  The  contour  plots  show  both  the 
frequency  domain  and  time  domain  of  the  results. 

First,  the  seven  carrier  frequencies  and  the  entire  cross  terms  are  shown  in  every 
figure.  Since  the  carrier  frequencies  are  1kHz,  2kHz,  3kHz,  4kHz,  5kHz,  6kHz  and 
7kHz,  the  cross  terms  will  be 

|{(l+2),  (1+3),  (1+4),  (1+5),  (1+6),  (1+7) 

,(2+3), (2+4), (2+5), (2+6), (2+7) 

,(3+4), (3+5), (3+6), (3+7) 

,(4+5), (4+6), (4+7) 

,(5+6), (5+7) 

,(6+7)}  kHz 

From  the  above  calculation,  the  cross  terms  appear  at  1.5  kHz  one  time,  2  kHz  one  time, 
2.5  kHz  two  times,  3  kHz  two  times,  3.5  kHz  three  times,  4  kHz  three  times,  4.5  kHz 
three  times,  5  kHz  two  times,  5.5  kHz  two  times,  6  kHz  one  time,  6.5  kHz  one  time.  In 
Figure  349  to  360,  only  the  1kHz  and  7kHz  are  real  carrier  frequencies,  all  the  others  are 
possibly  carrier  frequencies  or  cross  terms. 

Secondly,  in  Figure  56(b)  to  59(b)  there  marked  with  the  blue  arrows,  which 

indicate  the  frequency  sequence  from  the  top  to  the  bottom.  Take  Figure  56(b)  as  an 

example,  the  carrier  frequency  sequence  is,  as  the  blue  arrows  shown  from  the  top  to  the 

bottom,  AkIIz  — >  1kHz  — »  1  kHz  — >  6 kHz  — >  5 kHz  — >  2kHz  — >  2kHz . 
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As  can  be  seen,  the  colors  of  carrier  frequencies  and  cross  terms  in  Figure  57(b) 
are  different.  The  color  of  the  cross  terms  is  red  line  enclosed  by  yellow,  but  the  real 
carrier  frequencies  display  different  color  and  shapes.  They  are  yellow  circles  with  green 
curves  on  the  right  and  left  hand  sides.  According  to  Figure  10,  the  frequency  of  the 
cross  terms  is  half  of  the  two  real  carrier  frequencies.  The  time  period  of  the  cross  terms 
is  also  half  of  the  time  period  sum  of  the  two  real  carrier  frequencies. 

Taking  a  look  of  the  4kHz  in  Figure  57(b),  there  are  four  4kHz  bars  on  this  plot. 
From  the  top  to  the  bottom,  the  first  yellow  bar  is  the  real  carrier  frequency.  The  second 
red  bar  is  a  cross  tenn  caused  by  7kHz  and  1kHz,  the  appearing  time  period  is  right  on 
the  center  of  those  of  in  7kHz  and  1kHz.  The  third  red  bar  is  a  cross  tenn  caused  by 
6kHz  and  2kHz,  and  with  time  period  in  the  center  of  those  in  6kHz  and  2kHz.  The  forth 
red  bar  is  a  cross  tenn  caused  by  5kHz  and  3kHz  with  time  period  in  the  center  of  those 
in  5kHz  and  3kHz.  Note  there  are  two  3kHz  red  bars  in  Figure  57(b).  The  first  is  caused 
by  4kHz  and  2kHz  carrier  frequencies,  the  second  one  is  caused  by  1kHz  and  5kHz 
carrier  frequencies.  Another  3kHz  yellow  bar  is  a  real  carrier  frequency.  Figure  58  and 
Figure  59  show  another  sequence,  the  carrier  frequencies  are  also  examined  in  the 
figures. 


In  these  figures  the  time  scale  depends  on  the  sampling  frequency  and  the  number 


rf  sarrmle 

of  sample  points.  The  total  time  is  — - ,  which  can  be  used  to  determine  the  time  in 

f s 


each  frequency.  Since  the  total  sample  number  is  too  big  to  use  the  WD,  in  this  chapter 
WD  processes  350  points  for  Cycles/Frequency  =  10  and  700  points  for 
Cycles/Frequency  =  20. 
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(b) 

Figure  56.  WD  for  COSTAS  Frequency  Flopping,  Signal  Only,  (No. 85)  (a)  2D  Mesh  in 

Frequency  Domain,  (b)  Contour. 


110 


mesh  of  0^5^ 
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contour  of  C11520s 


(b) 

Figure  57.  WD  for  COSTAS  Frequency  Flopping,  Signal  Only,  (No. 88)  (a)  2D  Mesh  in 

Frequency  Domain,  (b)  Contour. 
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mesh  of  C21 71 0S 
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contour  of  C21710s 
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Figure  58.  WD  for  COSTAS  Frequency  Flopping,  Signal  Only,  (No. 91)  (a)  2D  Mesh  in 

Frequency  Domain,  (b)  Contour. 
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Figure  59.  WD  for  COSTAS  Frequency  Flopping,  Signal  Only,  (No. 94)  (a)  2D  Mesh  in 

Frequency  Domain,  (b)  Contour. 
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c. 


SUMMARY 


To  analyze  a  signal  using  Wigner  distribution,  the  cross  terms  are  unavoidable  for 
multi-frequencies  as  well  as  frequency-hopping  signals,  like  Costas.  The  cross  terms 
have  the  same  magnitude  as  the  real  carrier  frequencies.  The  frequency  of  any  cross  term 
is  half  of  the  sum  of  any  two  real  carrier  frequencies.  For  example,  if  and  fc2  are  real 


carrier  frequencies,  then  the  cross  term  generated  by  these  two  frequencies  is 


Both  real  carrier  frequencies  and  cross  term  exist  in  the  mesh  plots  and  contour 
plots.  Fortunately,  the  contour  plots  still  can  see  the  sequence  of  the  Costas  frequency 
hopping  signals.  As  long  as  the  carrier  frequencies  and  their  sequence  have  been  found, 
then  the  COSTAS  frequency  hopping  signals  can  be  recovered. 

Table  18  below  shows  the  detection  effectiveness  of  the  carrier  frequency 
sequence,  the  time  in  each  sequence  and  the  code  period  from  Figure  56  to  59.  They  are 
done  by  the  visible  inspection  of  the  WD  results. 
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COSTAS 

No. 

Sequence 

Time  in  Each  Frequency 

Code  Period 

85 

100% 

100% 

100% 

86 

80% 

50% 

50% 

87 

0% 

0% 

0% 

88 

100% 

100% 

100% 

89 

80% 

70% 

70% 

90 

0% 

0% 

0% 

91 

100% 

100% 

100% 

92 

80% 

50% 

50% 

93 

0% 

0% 

0% 

94 

100% 

100% 

100% 

95 

80% 

70% 

70% 

96 

0% 

0% 

0% 

Costas  code  effectiveness  on  detection 


Parameters 


□  Only  signal 

□  0  dB 

□  (-)  6  dB 


Table  18.  WD  Detection  Effectiveness  for  COSTAS  Code  Signals. 
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VII.  PHASE  SHIFT  KEYING/FREQUENCY  SHIFT  KEYING 

(PSK/FSK) 

A.  PSK/FSK  USING  A  COSTAS-BASED  FREQUENCY-HOPPING 

TECHNIQUE 

The  PSK/FSK  using  Costas  modulation  technique  is  the  result  of  a  combination  of 
frequency-shift  keying  based  on  a  Costas  frequency-hopping  matrix  and  phase-shift 
keying  using  Barker  sequences  with  different  lengths.  Here  will  briefly  describe  the 
phase  encoding  applied  to  a  Costas  signal,  generating  the  FSK/PSK  combined  waveform. 

Consider  a  Costas  frequency-hopped  signal,  the  firing  order  defines  which 
frequencies  will  appear  and  with  which  duration.  For  CW  radars,  the  usual  terminology 
does  not  apply  to  this  case.  Instead  of  a  “burst”  of  pulses,  frequencies  have  being 
continuously  emitted  during  a  defined  period  of  time.  This  period  may  be  divided  into 
sub-periods,  labeled  TF  for  each  frequency.  The  length  of  each  sub-period  depends  on 
the  sampling  interval.  During  each  sub-period,  as  the  signal  stays  in  one  of  the 
frequencies,  a  binary  phase  modulation  occurs  according  to  a  Barker  sequence  of  length 

five  (+  +  +  -  +),  seven  (+  +  +  --  +  -),  eleven  (+  +  + - +  --  +  -)  or  thirteen  (+  +  +  +  +  -- 

+  +  -  +  -  +).  [Ref.  23] 

For  example,  the  FSK/PSK  signal  defined  by:  S  =  {1+,  1+,  1+,  F,  1+,  2+,  2+,  2+,  2', 
2+,  3+,  3+,  3+,  3',  3+,  4+,  4+,  4+,  4',  4+,  5+,  5+,  5+,  5',  5+  }  represents  a  waveform  comprised 
of  five  different  frequencies,  which  are  each  subdivided  into  five  phase  slots,  labeled  Tp , 
according  to  the  Barker  sequence  of  length  five  (+  +  +  -  +). 
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The  final  waveform  may  be  seen  as  a  binary  phase-shifting  modulation  within 
each  frequency  hop,  resulting  in  25  phase  slots  equally  distributed  in  each  frequency. 

If  consider  NF  as  the  number  of  frequency  hops  and  Np  as  the  number  of  phase 
slots  of  duration  Tp  in  each  frequency  sub-period  TF ,  the  total  number  of  phase  slots  in 
the  FSK/PSK  waveform  is  given  by:  N  =  NF  ■  NP.  [Ref.  23] 


The  Barker  sequence  is  generated  and  the  frequency-hopping  signal  is  then  phase- 
modulated  accordingly.  For  example,  if  the  first  Costas  sequence  is  selected,  after  a 
phase  modulation  using  a  Barker  sequence  of  length  5,  the  final  waveform  becomes: 


S  =  4+,  4+,  4+,  4",  4+,  7+,  7+,  7+,  7",  7+,  1+,  1 1+,  1",  1+,  6+,  6+,  6+,  6",  6+,  5+,  5+,  5+, 


5',  5+,  2+,  2+,  2+,  2',  2+,  3+,  3+,  3+,  3',  3+  [Ref.  23] 


B.  FSK/  PSK  COMBINED  USING  A  TARGET-MATCHED  FREQUENCY 
HOPPING 


Instead  of  spreading  the  energy  of  the  signal  equally  over  a  broad  bandwidth,  this 
type  of  technique  concentrates  the  signal  energy  in  specific  spectral  locations  of  most 
importance  for  the  radar  and  its  typical  targets,  within  the  broad-spectrum  bandwidth. 
The  produced  signals  have  a  pulse  compression  characteristic,  and  therefore  they  can 
achieve  a  low  probability  of  intercept. 

The  implementation  starts  with  a  simulated-target  time-radar  response.  This  data 
is  the  Fourier  transformed,  the  frequency  components,  their  correspondent  magnitudes, 
and  their  initial  phases.  A  random  selection  process  chooses  each  frequency  with  a 
probability  distribution  function  defined  by  the  spectral  characteristics  of  the  target  of 

interest  obtained  from  the  Fast  Fourier  Transform  so  that  frequencies  at  the  spectral  peaks 
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of  the  target  (highest  magnitudes)  are  transmitted  more  often.  Each  “frequency  hop,” 
transmitted  during  a  specific  period  of  time,  is  also  modulated  in  phase,  having  its  initial 
phase  value  modified  by  a  pseudo-random  spreading-phase  sequence  code  of  values 
equally  likely  to  be  zero  or  n  radians. 

The  matched  FSK/PSK  radar  will  then  use  a  correlation  receiver  with  a  phase 
mismatched  reference  signal  instead  of  a  perfectly  phase  matched  reference.  This  allows 
the  radar  to  generate  signals  that  can  match  a  target’s  spectral  response  in  both  magnitude 
and  phase. 

This  simulation  serves  the  purpose  of  testing  the  perfonnance  of  a  new  digital 
cyclostationary  receiver  against  these  kinds  of  signals.  A  better  implementation  would 
require  an  iterative  solution  since  the  response  of  a  target  is  not  always  known 
beforehand.  [Ref.  23] 

C.  PSK/FSK  AND  TEST  SIGNALS  PAF  AND  PSD  ANALYSIS 

Eight  PSK/FSK  COSTAS,  eight  PSK/FSK  target  and  two  single  carrier  frequency 
signals  are  introduced  in  this  chapter.  They  are  numbered  from  97  to  114  and  listed  in 
Table  19.  The  “Sequence”  means  the  frequency  hopping  sequence.  The  numbers  in  the 
sequence  are  the  order  of  the  carrier  frequencies.  The  unit  of  these  carrier  frequencies  is 
kHz.  The  “Cycles/  Phase”  means  the  number  of  cycles  per  phase,  the  number  1,  5  and  10 
are  used  here.  The  sampling  frequency  is  15kHz  for  these  PSK/FSK  signals,  7kHz  for 
the  two  test  signals. 
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No. 

FSK/PSK  COSTAS 

Sequence 

Barker 

Code 

Length 

Cycles/ 

Phase 

SNR 

97 

FSK  PSK  C  1  15  5  5  s 

4716523 

5 

5 

Signal  Only 

98 

FSK  PSK  C  1  15  5  5  0 

4716523 

5 

5 

OdB 

99 

FSK  PSK  C  1  15  5  1  s 

4716523 

5 

1 

Signal  Only 

100 

FSK  PSK  C  1  15  5  1  0 

4716523 

5 

1 

OdB 

101 

FSK  PSK  C  1  15  11  5  s 

4716523 

11 

5 

Signal  Only 

102 

FSK  PSK  C  1  15  11  5  0 

4716523 

11 

5 

OdB 

103 

FSK  PSK  C  1  15  11  1  s 

4716523 

11 

1 

Signal  Only 

104 

FSK  PSK  C  1  15  11  1  0 

4716523 

11 

1 

OdB 

No. 

FSK/PSK  Target 

Sequence 

Random 

Hops 

Cycles/ 

Phase 

SNR 

105 

FSK  PSK  T  15  128  5  s 

1 

128 

5 

Signal  Only 

106 

FSK  PSK  T  15  128  5  0 

1 

128 

5 

OdB 

107 

FSK  PSK  T  15  256  5  s 

1 

256 

5 

Signal  Only 

108 

FSK  PSK  T  15  256  5  0 

1 

256 

5 

OdB 

109 

FSK  PSK  T  15  128  10  s 

2 

128 

10 

Signal  Only 

110 

FSK  PSK  T  15  128  10  0 

2 

128 

10 

OdB 

111 

FSK  PSK  T  15  256  10  s 

2 

256 

10 

Signal  Only 

112 

FSK  PSK  T  15  256  10  0 

2 

256 

10 

OdB 

No. 

Test  Signal 

fc  (kHz) 

fs  (kHz) 

SNR 

113 

T  1  7  1  s 

1 

7 

Signal  Only 

114 

T  2  7  1  s 

1  and  2 

7 

Signal  Only 

Table  19.  FSK/PSK  and  Test  Signals. 


D.  WIGNER  DISTRIBUTION  FOR  PSK/FSK  AND  TEST  SIGNALS 

Consider  the  same  eight  PSK/FSK  COSTAS,  eight  PSK/FSK  target  and  two 
single  carrier  frequency  signals  introduced  in  Table  19  and  use  them  as  the  inputs  for 
Wigner  Distribution  (WD).  The  WD  results  of  these  twelve  signals  are  showed  in  Figure 
60  to  67.  The  mesh  plots  show  the  frequency  domain  of  these  signals  after  WD.  The 
contour  plots  show  both  the  frequency  domain  and  time  domain  of  the  results. 

The  following  figures  for  the  PSK/FSK  COSTAS  signals  are  very  similar  to  the 


figures  of  COSTAS  signals.  Comparing  with  Figure  60(b)  and  Figure  56(b)  and  57(b), 
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the  Wigner  distribution  analysis  gets  very  similar  results  for  the  PSK/FSK  COSTAS 
signals  and  COSTAS  signals.  The  carrier  frequencies  sequence,  4716523,  can  be  seen  as 
yellow  circles  and  the  frequency  cross  terms  can  be  seen  as  red  bars  in  Figure  60(b).  As 
a  matter  of  fact,  no  matter  whether  the  5  bits  or  1 1  bits  Barker  code  length,  the  time- 
frequency  relationships  are  the  same  for  the  PSK/FSK  COSTAS  signals. 

In  Figure  61(b)  and  63(b),  the  CPP  is  one,  so  the  shape  of  the  carrier  frequencies 
are  different  from  those  with  five  cycles  per  phase  in  Figure  60(b)  and  62(b).  One  can 
still  distinguish  the  real  carrier  frequencies  by  the  lighter  colors. 

Since  the  phase  changes  are  by  random,  for  the  PSK/FSK  target  signals  neither 
the  carrier  frequencies  nor  the  cross  terms  can  be  examined  on  the  mesh  plots  and 
contour  plots.  Figure  68  shows  the  one  carrier  frequency  test  signal.  Figure  69  shows  the 
two  carrier  frequencies  and  frequency  cross  term  as  described  in  Chapter  II. 
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Figure  60.  WD  for  FSK/PSK  COSTAS  with  5  bit  Barker  Code,  CPP  =  5,  Signal  Only, 
(No. 97)  (a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  61.  WD  for  FSK/PSK  COSTAS  with  5  bit  Barker  Code,  CPP  =  1,  Signal  Only, 
(No. 99)  (a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  62.  WD  for  FSK/PSK  COSTAS  with  1 1  bit  Barker  Code,  CPP  =  5,  Signal  Only, 
(No.  101)  (a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  63.  WD  for  FSK/PSK  COSTAS  with  1 1  bit  Barker  Code,  CPP  =  1,  Signal  Only, 
(No.  103)  (a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  64.  WD  for  FSK/PSK  Target  with  128  Random  Hops,  CPP  =  5,  Signal  Only, 
(No.  105)  (a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  65.  WD  for  FSK/PSK  Target  with  256  Random  Hops,  CPP  =  5,  Signal  Only, 
(No.  107)  (a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  66.  WD  for  FSK/PSK  Target  with  128  Random  Hops,  CPP  =  10,  Signal  Only, 
(No.  109)  (a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  67.  WD  for  FSK/PSK  Target  with  256  Random  Hops,  CPP  =  10,  Signal  Only, 
(No. Ill)  (a)  2D  Mesh  in  Frequency  Domain,  (b)  Contour. 
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Figure  68.  WD  for  1  kHz  Carrier  Frequency  Test  Signal,  (No.l  13)  (a)  2D  Mesh  in 

Frequency  Domain,  (b)  Contour. 
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Mesh  of  two  carrier  frequencies  example 
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Contour  of  two  carrier  frequencies  example 
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Figure  69.  WD  for  1kHz  and  2  kHz  Carrier  Frequency  Test  Signal,  (No.  1 14)  (a)  2D 

Mesh  in  Frequency  Domain,  (b)  Contour. 
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E.  SUMMARY 

From  the  above  Figures,  only  the  PSK/FSK  COSTAS  signals  and  the  last  two  test 
signals  show  the  carrier  frequencies  or  their  sequence.  Wigner  Distribution  cannot 
analyze  the  PSK/FSK  target  frequency  hopping  signals. 

Table  20  below  shows  the  detection  effectiveness  of  the  sequence,  bandwidth, 
code  period  and  phase  from  Figure  60  to  67. 


FSK/PSK_COSTAS 

No. 

Sequence 

Bandwidth 

Code  Period 

97 

100% 

100% 

100% 

98 

50% 

50% 

50% 

99 

100% 

100% 

100% 

100 

50% 

50% 

50% 

101 

100% 

100% 

100% 

102 

50% 

50% 

50% 

103 

100% 

100% 

100% 

104 

50% 

50% 

50% 

FSK/PSKTARGET 

No. 

Sequence 

Bandwidth 

Phase 

105 

0% 

0% 

0% 

106 

0% 

0% 

0% 

107 

0% 

0% 

0% 

108 

0% 

0% 

0% 

109 

0% 

0% 

0% 

110 

0% 

0% 

0% 

111 

0% 

0% 

0% 

112 

0% 

0% 

0% 
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FSK/PSK  Costas  effectiveness  on  detection 


□  Only  signal 
■  0  dB 


Parameters 


FSK/PSK  Target  effectiveness  on  detection 


Parameters 


□  Only  signal 
■  0  dB 


Table  20.  WD  Detection  Effectiveness  for  FSK/PSK  Signals. 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Nine  types  of  Low  Probability  of  Intercept  (LPI)  signals  were  considered  in  this 
thesis.  They  are  Binary  Phase  Shift  Keying  (BPSK),  Frequency  Modulation  Continuous 
Wave  (FMCW),  Frank  code,  PI  code,  P2  code,  P3  code,  P4  code,  COSTAS  frequency 
hopping  and  Phase  Shift  Keying/Frequency  Shift  Keying  (PSK/FSK)  signals.  There  are 
114  signals  were  generated  for  Wigner  distribution  analysis  in  this  thesis.  The  detail 
parameters  of  these  114  test  signals  are  listed  in  Appendix  B. 

The  Wigner  distribution  is  a  two-dimension  function  describing  the  frequency 
content  of  a  signal  as  a  function  of  time.  The  Wigner  distribution  has  been  noted  as  one 
of  the  more  useful  time-frequency  analysis  techniques  for  signal  processing.  Using  the 
Wigner  distribution  one  can  identify  the  frequency  and  time  changes  in  most  of  the  LPI 
signals.  This  thesis  successfully  examined  the  nine  types  of  LPI  signals  mentioned 
above.  All  the  signal  parameters  excluding  the  PSK/FSK  signals  can  be  represented  by 
Wigner  distribution.  For  those  detectable  signals,  some  of  them  can  show  the  parameters 
at  both  SNR  =  OdB  and  -6  dB,  some  of  them  cannot  be  detected  at  SNR  =  -6dB.  That  is 
to  say  the  signal  to  noise  ratio  is  a  main  concern  to  the  analysis  of  the  LPI  signals.  Please 
see  the  technical  report.  [Ref.  7] 

Other  than  the  signal  to  noise  ratio,  there  are  some  other  important  characteristics 
that  need  to  be  emphasized.  For  BPSK  signals,  the  number  of  Barker  bits  can  be 
displayed  clearly  only  when  the  number  of  cycles  per  barker  bit  is  one  ( NPBB  =  1).  For 
FMCW  signals,  the  parameters  can  be  displayed  meaningfully  when  the  modulation 
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period  less  than  or  equal  to  half  of  the  sample  points  divided  by  sampling  frequency 


(X  ^ 


#  sample 

2 7, 


)  and  when  the  modulation  bandwidth  less  than  or  equal  to  2  times  of 


carrier  frequency  (  AF  <  2fc )  or  (  AF  <  2{fs  /2-fc)  ).  For  the  poly  phase  signals 
including  Frank  codes,  PI  codes,  P2  codes,  P3  codes  and  P4  codes,  the  phase  length  (N or 
N2)  can  be  obtained  precisely  when  the  number  of  cycles  per  phase  is  one  ( CPP  =  1). 
The  higher  the  CPP ,  the  harder  it  is  to  get  the  right  phase  length. 


Last,  the  frequency  cross  terms  are  irremovable  for  COSTAS  frequency  hopping 
signals  or  any  other  multi-carrier-frequency  signal.  The  cross  terms  have  the  same 
amplitude  as  real  carrier  frequencies  and  their  frequencies  are  the  average  of  any  two  real 

carrier  frequencies  ( ^cl  cl  ). 

B.  RECOMMENDATIONS 


The  Choi-Williams  Distribution  (CWD)  is  a  time  frequency  analysis  technique 
that  evolved  from  Wigner  Distribution  (WD).  The  most  important  advantage  of  using  the 
Choi-Williams  distribution  is  that  it  has  the  ability  to  eliminate  the  frequency  cross  tenns. 
The  definition  equation  of  the  Choi-Williams  Distribution  was  introduced  in  chapter  I, 
but  no  experiments  were  examined  in  this  thesis.  The  CWD  is  deserving  of  further  study. 

The  output  of  the  Wigner  distribution  is  a  real  two-dimensional  matrix  with  the 
size  of  length  (/)  times  length  (0,  where  /  and  Q  are  the  in-phase  and  the  quadrature 
component  of  the  input  signal.  The  input  signal  in  this  thesis  was  defined  as  x  =  I  +  j*Q. 
If  the  length  of  /  or  Q  is  10k  points,  then  the  size  of  the  output  matrix  will  be  100M 
points.  In  other  words,  the  output  matrix  will  occupy  800MB  memory. 
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Thus,  the  bigger  the  input  signal  size,  the  better  and  more  information  will  be 
shown  in  the  mesh  and  contour  plots.  But  the  size  of  the  output  matrix  is  too  great  for 
MATLAB®  processing  on  the  computer  used  in  this  study;  MATLAB®  ran  out  of 

memory. 

The  specifications  of  the  computer  used  to  run  the  Wigner  distribution 
MATLAB®  program  is  Intel  Pentium©  4  CPU  with  2.0  GHz  clock  speed  and  1  GB 
RAM.  It  is  representative  of  the  best  in  the  personal  computer  market.  Future  Wigner 
distribution  analysis  obviously  requires  more  capable  computers. 
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APPENDIX  A.  MATLAB  CODE 


A.  WIGNERM.M 


function  [W]  =  wigner_fft (x) 

%  Design  a  variable  W  as  the  Wigner  distribution  function  for 
input  signal  x 

%  Paul  Milne 

%  EC4700  Project  in  Naval  Postgraduate  School 
%  May  2001 
"6 

%  Modified  by  Jen-Yu,  Gau 
%  Jan.  2002 

N  =  length (x)/2; 
input  signal  x 
odd  =  0 ; 
even 

if  mod(2*N,2)  ~=  0  %  If  the  signal  has  odd  length 

disp (' WARNING :  Input  Sequence  Odd:'); 

disp  ( '  A  zero  will  be  appended  for  calculations.'); 

%  Display  warning  messages 

x  =  [0  x] ;  %  Add  a  zero  at  the  beginning  to  make  even  length 

N  =  length (x)/2;  %  Re-compute  length  2N 

odd  =1;  %  Mark  the  original  input  signal  has  odd  length 

end 

M  =  2*N-1 
a  =  -N; 
b  =  N-l ; 
y  =  x; 
signal 

for  n  =  a:b  %  Set  index  n  from  -N  to  N-l 

x  =  shiftz(y,n);  %  Call  function  shiftz.m  to  shift  the 

initial  input  signal  x 

xl  =  [x ( 1 , N :  M)  0  x ( 1 , 1 : N-l ) ] ; 

x2  =  [fliplr (x (1, 1 :N) )  0  f liplr (x ( 1 , N+l :  M) ) ] ; 

%  Compute  the  shifted  signal  xl  and  x2  in  order  to  get  the 
kernel 

kern  =  xl . *con j  (x2 ) ;  %  Compute  the  kernel  as  in  Eq(2.5) 

W(n-a+l,:)  =  fft (kern) ; 

%  Apply  the  fast  Fourier  transform  (FFT)  algorithm  of  Matlab 
command 

end 

W=2*real (W)  %Wigner  distribution  results  are  always  real 

%  The  final  result  of  the  Wigner  distribution  as  in  Eq(2.13) 

if  odd  %  If  the  original  input  signal  has  odd  length 


%  Set  M  =  2N-1 

%  Set  a  as  integration  down  limit 
%  Set  b  as  integration  up  limit 

%  Set  y  as  another  name  for  the  input 


%  Set  2N  as  the  length  of  the 
%  Initial  the  signal  length  as 
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%  Then  remove  extra  output  created  by  the  adding  zero 
W  =  W (2 : 2*N, 1 : 2*N-1 ) ;  %  The  final  result  becomes 

end 


B.  SHIFTZ.M 


function  [y]  =  shift (x,n) 

%  The  function  is  used  to  shift  input  signal  x  with  n  zeros  to  the 
right  and  left 

%  in  order  to  get  the  kernel  as  in  Eq(2.5) 

%  When  x  length  2N,  then  n  is  from  -N  to  N-l 

%  Paul  Milne 

%  EC4700  Naval  Postgraduate  School 
%  May  2001 
"6 

%  Modified  by  Jen-Yu,  Gau 
%  Jan.  2002 

N  =  length (x)/2;  %  Set  2N  as  the  length  of  the  input  signal  x 

if  n  >=  0  %  First  shift  the  up  half  part  of  x 

y  =  [zeros (l,n)  x(l:2*N-n)];  %  Shift  x  to  the  right  with  n 

zeros 

else  %  Then  shift  the  down  half  part  of  x 

n  =  abs (n) ; 

y  =  [x(l+n:2*N)  zeros (1, n) ] ; 

%  Shift  x  to  the  left  with  n  zeros 
end 
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APPENDIX  B.  TEST  MATRIX 


No. 

BPSK 

Barker  Code 
Length 

NPBB 

SNR 

1 

B  1  7  7  1  s 

7 

1 

Signal  Only 

2 

B  1  7  7  1  0 

7 

1 

OdB 

3 

B  1  7  7  1  -6 

7 

1 

-6  dB 

4 

B  1  7  11  1  s 

11 

1 

Signal  Only 

5 

B  1  7  11  1  0 

11 

1 

OdB 

6 

B  1  7  11  1  -6 

11 

1 

-6  dB 

7 

B  1  7  7  5  s 

7 

5 

Signal  Only 

8 

B  1  7  7  5  0 

7 

5 

OdB 

9 

B  1  7  7  5  -6 

7 

5 

-6  dB 

10 

B  1  7  11  5  s 

11 

5 

Signal  Only 

11 

B  1  7  11  5  0 

11 

5 

OdB 

12 

B  1  7  11  5  -6 

11 

5 

-6  dB 

FMCW 

Mod.  BW 
(Hz) 

Mod.  Period 
(ms) 

SNR 

13 

F  1  7  250  20  s 

250 

20 

Signal  Only 

14 

F  1  7  250  20  0 

250 

20 

OdB 

15 

F  1  7  250  20  -6 

250 

20 

-6  dB 

16 

F  1  7  250  30  s 

250 

30 

Signal  Only 

17 

F  1  7  250  30  0 

250 

30 

OdB 

18 

F  1  7  250  30  -6 

250 

30 

-6  dB 

19 

F  1  7  500  20  s 

500 

20 

Signal  Only 

20 

F  1  7  500  20  0 

500 

20 

OdB 

21 

F  1  7  500  20  -6 

500 

20 

-6  dB 

22 

F  1  7  500  30  s 

500 

30 

Signal  Only 

23 

F  1  7  500  30  0 

500 

30 

OdB 

24 

F  1  7  500  30  -6 

500 

30 

-6  dB 

FRANK 

N  (Phase) 

Cycles/Phase 

SNR 

25 

FR  1  7  4  1  s 

4 

1 

Signal  Only 

26 

FR  1  7  4  1  0 

4 

1 

OdB 

27 

FR  1  7  4  1  -6 

4 

1 

-6  dB 

28 

FR  1  7  4  5  s 

4 

5 

Signal  Only 

29 

FR  1  7  4  5  0 

4 

5 

OdB 

30 

FR  1  7  4  5  -6 

4 

5 

-6  dB 

31 

FR  1  7  8  1  s 

8 

1 

Signal  Only 

32 

FR  1  7  8  1  0 

8 

1 

OdB 

33 

FR  1  7  8  1  -6 

8 

1 

-6  dB 

34 

FR  1  7  8  5  s 

8 

5 

Signal  Only 

35 

FR  1  7  8  5  0 

8 

5 

OdB 

36 

FR  1  7  8  5  -6 

8 

5 

-6  dB 

PI 

N  (Phase) 

Cycles/Phase 

SNR 
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37  PI  1  7  4  1  s 

38  PI  1  7  4  1  0 

39  PI  1  7  4  1  -6 

40  PI  1  7  4  5  s 

41  PI  1  7  4  5  0 

42  PI  1  7  4  5  -6 

43  PI  1  7  8  1  s 

44  PI  1  7  8  1  0 

45  PI  1  7  8  1  -6 

46  PI  1  7  8  5  s 

47  PI  1  7  8  5  0 

48  PI  1  7  8  5  -6 

_ P2 

49  P2  1  7  4  1  s 

50  P2  1  7  4  1  0 

51  P2  1  7  4  1  -6 

52  P2  1  7  4  5  s 

53  P2  1  7  4  5  0 

54  P2  1  7  4  5  -6 

55  P2  1  7  8  1  s 

56  P2  1  7  8  1  0 

57  P2  1  7  8  1  -6 

58  P2  1  7  8  5  s 

59  P2  1  7  8  5  0 

60  P2  1  7  8  5  -6 

_ P3 

61  P3  1  7  16  1  s 

62  P3  1  7  16  1  0 

63  P3  1  7  16  1  -6 

64  P3  1  7  16  5  s 

65  P3  1  7  16  5  0 

66  P3  1  7  16  5  -6 

67  P3  1  7  64  1  s 

68  P3  1  7  64  1  0 

69  P3  1  7  64  1  -6 

70  P3  1  7  64  5  s 

71  P3  1  7  64  5  0 

72  P3  1  7  64  5  -6 

_ P4 

73  P4  1  7  16  1  s 

74  P4  1  7  16  1  0 

75  P4  1  7  16  1  -6 

76  P4  1  7  16  5  s 

77  P4  1  7  16  5  0 

78  P4  1  7  16  5  -6 


N  (Phase) 
4 
4 
4 
4 
4 
4 
8 
8 
8 
8 
8 
8 

N  (Phase) 
16 
16 
16 
16 
16 
16 
64 
64 
64 
64 
64 
64 

N  (Phase) 
16 
16 
16 
16 
16 
16 


Cycles/Phase 

1 

1 

1 

5 

5 

5 

1 

1 

1 

5 

5 

5 

Cycles/Phase 

1 

1 

1 

5 

5 

5 

1 

1 

1 

5 

5 

5 

Cycles/Phase 

1 

1 

1 

5 

5 

5 


Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 
SNR 

Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 
SNR 

Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 
SNR 

Signal  Only 
OdB 
-6  dB 

Signal  Only 
OdB 
-6  dB 


79 

P4  1  7  64  1  s 

64 

1 

Signal  Only 

80 

P4  1  7  64  1  0 

64 

1 

OdB 

81 

P4  1  7  64  1  -6 

64 

1 

-6  dB 

82 

P4  1  7  64  5  s 

64 

5 

Signal  Only 

83 

P4  1  7  64  5  0 

64 

5 

OdB 

84 

P4  1  7  64  5  -6 

64 

5 

-6  dB 

COSTAS 

Sequence 

Cycles/Freq. 

SNR 

85 

C  1  15  10  s 

4716523 

10 

Signal  Only 

86 

C  1  15  10  0 

4716523 

10 

OdB 

87 

C  1  15  10  -6 

4716523 

10 

-6  dB 

88 

C  1  15  20  s 

4716523 

20 

Signal  Only 

89 

C  1  15  20  0 

4716523 

20 

OdB 

90 

C  1  15  20  -6 

4716523 

20 

-6  dB 

91 

C  2  17  10  s 

2638751 

10 

Signal  Only 

92 

C  2  17  10  0 

2638751 

10 

OdB 

93 

C  2  17  10  -6 

2638751 

10 

-6  dB 

94 

C  2  17  20  s 

2638751 

20 

Signal  Only 

95 

C  2  17  20  0 

2638751 

20 

OdB 

96 

C  2  17  20  -6 

2638751 

20 

-6  dB 

FSK/PSK  COSTAS 

Sequence 

Barker 

Code 

Length 

Cycles/ 

Phase 

SNR 

97 

FSK  PSK  C  1  15  5  5  s 

4716523 

5 

5 

Signal  Only 

98 

FSK  PSK  C  1  15  5  5  0 

4716523 

5 

5 

OdB 

99 

FSK  PSK  C  1  15  5  1  s 

4716523 

5 

1 

Signal  Only 

100 

FSK  PSK  C  1  15  5  1  0 

4716523 

5 

1 

OdB 

101 

FSK  PSK  C  1  15  11  5  s 

4716523 

11 

5 

Signal  Only 

102 

FSK  PSK  C  1  15  11  5  0 

4716523 

11 

5 

OdB 

103 

FSK  PSK  C  1  15  11  1  s 

4716523 

11 

1 

Signal  Only 

104 

FSK  PSK  C  1  15  11  1  0 

4716523 

11 

1 

OdB 

FSK/PSK  Target 

Sequence 

Random 

Hops 

Cycles/ 

Phase 

SNR 

105 

FSK  PSK  T  15  128  5  s 

1 

128 

5 

Signal  Only 

106 

FSK  PSK  T  15  128  5  0 

1 

128 

5 

OdB 

107 

FSK  PSK  T  15  256  5  s 

1 

256 

5 

Signal  Only 

108 

FSK  PSK  T  15  256  5  0 

1 

256 

5 

OdB 

109 

FSK  PSK  T  15  128  10  s 

2 

128 

10 

Signal  Only 

110 

FSK  PSK  T  15  128  10  0 

2 

128 

10 

OdB 

111 

FSK  PSK  T  15  256  10  s 

2 

256 

10 

Signal  Only 

112 

FSK  PSK  T  15  256  10  0 

2 

256 

10 

OdB 

Test  Signal 

fc  (kHz) 

fs  (kHz) 

SNR 

113 

T  1  7  1  s 

1 

7 

Signal  Only 

114 

T  2  7  1  s 

1  and  2 

7 

Signal  Only 
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